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International Conference on Shock Waves in 
Condensed Matter 


957A1193A Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, p3 


[Article by Professor A.A. Borisov, chairman, confer- 
ence organizing Committee] 


{(FBIS Translated Text) Materials from an international 
conference on shock waves in condensed matter will 
be published in the January and February issues of this 
journal. The conference was held in St. Petersburg from 
18 through 22 July 1994. It covered all aspects of the 
generation, propagation, and interaction of shock waves 
in solids and liquids and phase and chemical transfor- 
mations occurring under high dynamic loads, including 
detonation processes and the transition from combus- 
lion to detonation, pressure wave-induced fracture and 
deformation, and characteristic equations. A total of 47 
oral reports and 105 podium reports were presented. A 
total of 185 individuals participated in the conference, 
including 148 scientists from CIS countries who repre- 
sented all of the leading scientific research centers at 
which research in the field is being conducted and 37 
scientists from France, the United States, England, and 
Germany. 


The conference's success was made possible by finan- 
cial support from the Russian Federation Ministry of 
the Nuclear Industry and Russian Federation Ministry 
of Science and Technical Policy. 


It was decided that the next conference would be held 
in the city of Samara in 2 years’ time. 


Flame Bench Tests To Determine Whether a 
Solid-Fuel-Fired Rocket Engine Can Withstand a 
Short-Term Lateral Load 


957A1201A Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 11-17 


[Article by A.V. Ostrik and V.P. Petrovskiy, Central 
Physicotechnical Institute, Russian Federation Ministry 
of Defense; manuscript received 20 Sep 94; first para- 
graph is KHIMICHESKAYA FIZIKA abstract; UDC 
§32.598:620.17:678.01] 


(FBIS Translated Text] An experimental method of 
Studying the consequences of brief intense loads on 
working solid-fuel-fired rocket engines is described. The 
results are presented from tests of four of the same 
type of model engine that were, during the course of 
their operation, subjected to lateral loading by pulse 
pressures with different amplitudes and durations. It 
is demonstrated that a nondesign increase in internal 
pressure as a result of an increase in combustion surface 
in cases of spalling of the fuel charge as a result of a 


pressure wave is one of the main reasons for fracture of 
the combustion chamber. 


Both when operated on the ground and in flight, a 
solid-fuel-fired rocket engine is subjected not just to 
Static loads but also to numerous nonstauonary loads 
of a shock and compressive nature that develop as 
a result of unevenness of the road along which they 
are being transported, pressure pulses and changes 
in internal pressure when the engine is started and 
stopped, wind gusts, the shock wave and fragmentat n 
generated during an explosion, and intense directed 
energy fluxes.'* Depending on their source, the duration 
of such loads ranges from 10° to 10° seconds, and their 
characteristic pressure pulses range from 0.1 kPas to 
5 kPas. An engine cannot be made to withstand loads 
with the aforesaid parameters without a detailed study of 
the processes occurring in an operating solid-fuel-fired 
rocket engine at the moment of and after application of 
the load. 


Nonstationary loading of a structure causes stratifica- 
tion, cracks, and other defects in the composite mater- 
als of the combustion chamber and nozzle block and, in 
cases of high pressure amplitudes, initiates breakdown 
and explosive processes in fuel. Modeling methods* can- 
not be used to study the development of fractures in bod- 
ies or explosive processes in fuel because the require- 
ments that the criterial parameters for the model and an 
actual engine virtually coincide with one another means 
that they must be identical from the standpoints of ab- 
solute dimensions and the properties of their materials 
It is, for example, altogether impossible to reproduce 
shock wave-initiated detonation under conditions where 
a fuel charge is geometrically similar to a model but 
has a diameter not exceeding the critical diameter of 
detonation even though such detonation can in princi- 
ple occur in an actual engine. Analogous problems also 
arise in attempts to model the fracture of a body based 
on a model. The condition of equality of the relative 
thicknesses A/R (where h and R are the thickness and 
characteristic radius of the shell) of a composite body 
and its model with the thickness of the reinforcing fil- 
ament remaining constant results in a decrease in the 
number of reinforcing layers in the model in relation to 
the actual body, which in turn distorts the nature and 
sequence of fracture of the said layers. 


A large number of theoretical’'’ and experimental'** 
Studies directly related to the dynamic strength of solid- 
fuel-fired rocket engines have been conducted to date 
Three of them’* have examined the effects of axisym- 
metric external pressure on a long cylindrical shel! with 
an isotropic or orthotropic filler. The propagation and re- 
flection of waves in the filler where its inner boundary 
is either firmly attached to a deformable body or free of 
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stress have been considered. In two studies'®'' the dis- 
tribuuon of the pressure pulse was also assumed to be 
axisymmetric, however, shells with a finite length were 
studied. One study'* examined infinite shells with a filler 
under the effect of a load moving along the generatrix. 
Another study'’ presented a method for calculating the 
effect of a nonaxisymmetric load of low duration on an 
infinitely long cylindrical shell with an elastic filler. 


The experimental studies have primarily been devoted 
to different aspects of the dynamic strength of solid- 
fuel-fired rocket engine components. One study’* has 
described experiments designed to determine the frac- 
ture of a solid-fuel-fired rocket engine with a firmly at- 
tached charge under the effect of an evenly distributed 
compressive load or internal pressure. The shock wave 
sensitivity of a mixed two-base fuel has been exam- 
ined,'* and data regarding the levels of spalling of solid 
rocket fuel have been obtained.'® Two studies!’ '* have 
been devoted to the characteristic features of the spalling 
of solid fuel charges with a complex inner channel ge- 
ometry during nonaxisymmetric external pulse loading. 


The studies that have examined the effect of nonstation- 
ary loads on the components of a nonoperating solid- 
fuel-fired rocket engine do not make it possible to deter- 
mine the safe levels of such loads because they provide 
virtualiy no data regarding the maximum dimensions 
and nature of damage incurred by each component be- 
fore the operation of the engine as a whole its inter- 
rupted. Another important question that remains open is 
the question of the effect that defects formed in a charge 
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after application of a load have on the rate of fuel com 
bustion and, consequently, on the increase in internal 
pressure that may result in fracture of the body's com- 
posite shell. Furthermore, the main causes of fracture 
of an engine during nonstationary lateral loading have 
not been definitively identificd either. To answer these 
questions it is necessary to conduct flame bench tests in 
which the engine is subjected to pressure pulses while 
it is Operating. Such tests make it possible to obtain 
detailed information about the functioning of a solid- 
fuel-fired rocket engine after it has been damaged by a 
nonstationary load. It is important that the engine's com- 
ponents be damaged under conditions that are as close 
to flight conditions as possible inasmuch as damage de- 
formation of the structure occurs when actual fuel (for 
Safety's sake, an inert filler is generally used in dynam 
tests) and internal pressure are present and the temper- 
ature of the fuel charge’s inner layers 1s elevated 


1. Explosive Devices for Creating Loads 


In the present stage, the gas dynamic method of local! 
loading has gained widespread popularity. The method 
entails creating a nonstationary load by a shock wave 
generated by detonating an explosive that is distributed 
three-dimensionally along the surface of the structure 
being tested or along a porous (rubber) spacer or on an 
inert substrate located a specified distance away.’ Table 
1 presents data about explosive devices for creating 
nonstationary loads that have been used to test operating 
solid-fuel-fired rocket engines. 





Table 1. Explosive Devices for Creating Nonstationary Loads 














Parameters 
Device Duration (7), jus Pressure Pulse (I), 
kPas 
Contact sector charge without spacer 12 ORS 
Contact sector charge with spacer made of porous rubber 5-10 0.1.2 





Equidistant surface charge 











10-100 04.3 





A contact sector charge 1s designed to create loads of 
microsecond duration. The charge is made from a sheet 
explosive that has been broken up into small segments/ 
sectors (with a gap of 0.5 mm _ between adjacent 
segments) to eliminate any knife effect’ (the said effect 
arises during the interaction of shock waves arriving 
from different charge initiation points) and is placed 
directly on the barrier or on a porous layer of rubber 


(p = 0.7 g/cm’) attached to the barner (the rubber, 
explosive, and barrier are glued together). The rubber 
layer makes it possible to increase the duration of 
loading to as much as 10 ys depending on the thickness 
of the spacer. Each sector of the charge is initiated at 
one point by the minimum strip of explosive possible 
for a stable detonation with a width of 0.5 cm, which 
virtually eliminates any local effect of the explosive at 
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the sites where the charge and strip have been glued. 
Strips of identical length are gathered into a bundle 
from one side and initiated by a high-speed detonator. 
The detonator is removed from the loading surface 
and does not create any load on the test objects that 
had not been taken into consideration when the charge 
was designed. Simultaneity of loading is achieved by 
mulupoint initiation. The law governing distribuuon 
of the pressure pulse along the solid-fuel-fired rocket 
engine's surface is realized by varying the thickness of 
the explosive or by perforating it. 


An equidistant surface charge is an elastic sheet explo- 
sive On an inert substrate that is located a specified dis- 
tance away from the test object (at a distance of 30 cm, 
the duration of loading equals 100 ys). The inert sub- 
Strate and charge are attached to a metal wire frame that 
rests on the structure undergoing testing. Simultaneity 
of detonation of the explosive is achieved by multipoint 
initiation analogous to that described for contact sector 
charges. 


After the calculation portion of the design process 
(i.e., after the distribution of the explosive’s mass in 
space) has been completed, freely suspended pendulum 
targets (whose shape and dimensions are similar to 
those of the test object) are used to work out the 
charges by experimentation. The total pulse created 
by the explosive device can be determined based on 
the angle of deviation of the pendulum targets. When 
an equidistant surface charge is being worked out, the 
distribution of pressure along the target’s surface and its 
change over ume are measured by piezoelectric sensors 
in addition to the total pulse. 


2. Measured Parameters 


Two groups of parameters were measured during study 
of the nonstationary processes occurring in solid-fuel- 
fired rocket engine operating under conditions of short- 
term lateral loads. The first group characterizes the 
engine's response during application of the load and 
includes nonstationary deformation of the body and 
overloads in important subassemblies of the solid-fuel- 
fired rocket engine. The second group of parameters 
characterizes the object's condition after it has been 
subjected to the load. The following parameters were 
selected as key parameters of the second group: internal 
pressure in the combustion chamber, the engine's thrust, 
and quasi-static deformations and temperature of the 
body. 


Nonstationary deformations were measured by elec- 
trotensometry by using a balanced single bridge. Type 
KB-1(0-200 strain gauge sensors were used as measuring 
sensors. They made it possible to measure relative de- 


formation to «. approximately equal to 4 percent with 
a Maximum error not exceeding 15 percent. 


A type ADP-10-01 piezoclectric converter with a high 
conversion coefficient was used to measure overloads 
It is a high-frequency, vibration-resistant converter that 
makes it possible to measure accelerations up to 10% pg 


3. Test Results 


The main difficulues in conducung strength tests of op- 
erating solid-fuel-fired rocket engines have to do with 
vibrations, heat flows from the jet, and the presence 
of a rigid jig that distorts the parameters characterizing 
the engine's response to a nonstationary load. Experi- 
ments have demonstrated that heat flows and vibrations 
do not significantly affect the error of the measuring 
instruments or the serviceability of the explosive de- 
vices. The jig was redesigned to reduce its rigidity: The 
thick frame in the center part of the jig was replaced 
in the 120° sector by a flexible metal tension strip. Aft- 
terward, the unfilled solid-fuel-fired rocket engine body 
was tested in the jig and without it (with the body freely 
suspended by a cable that was attached to the anterior 
bottomplate). The fracture pattern and process of the de- 
velopment and amplitude of nonstauonary deformations 
turned out to be similar in both cases (specifically, the 
axial deformations in the center of the region incurring 
the loading did not differ by more than 10 percent) 


First, the operating solid-fuel-fired rocket engine was 
subjected to the load of an equidistant surface charge 
with |, = 0.3 kPas and 7,= 100 jus, which judging 
by the results of preliminary numerical calculations,’ 
should not cause it to fracture. When the load was 
applied approximately to the center of the engine's 
operation, it had virtually no effect. Only slight changes 
in internal pressure and thrust were observed. It should 
be noted that in the center of the loading, nonstationary 
longitudinal tensile deformation was added to static 
deformation “stat approximately equal to 1.3 percent) 
and reached a level of ¢ approximately equal to 1.6 
percent, which 1s near the critical level. Significant 
overloads on the combustion chamber body and at 
the site where the nozzle block is attached were also 
observed while the load was being applied. A spectral 
analysis of the process of development of the overloads 
demonstrated that the vibration spectrum of the system 
body-charge has a distinct peak at a frequency of 600 
Hz and weaker peaks at frequencies of 1.6 and 2.6 kHz 


An engine undergoing testing thus possesses sufficient 
strength to withstand a load with I = ().3 kPas and 1, 
= 100 ys. At the same time, applying a pressure pulse 
to a solid-fuel-fired rocket engine from a contact sector 
charge through rubber with values of / = 2 kPas and 7, = 
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10 ys leads to fracture of the internal pressure-damaged 
body after a period somewhere on the order of 900 ys 
after initiation of the charge. At the moment when the 
body loses its load-bearing capability, the longitudinal 
deformations reach a level of € = 2 percent. No increase 
in intesnal pressure after application of the load was 
observed, and the main reason for disruption in the 
engine’s operauion was the decrease in the combustion 
chamber’s strength that occurred as a result of the pulse 
loading. 
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Further tests were conducted at lower pulse pressure 
levels (J, = 0.5 kPas and 7, = 10 yes) to detect causes of 
fracture of the solid-fuel-fired rocket engine other than 
damage to the body; however, subjecung the engine to 
a lesser pulse at the end of its operation did not lead to 
loss of the body’s seal without a preliminary increase 
in pressure either. 





Table 2. Pressure (in MPa) Inside the Combustion Chamber of an Operating Engine 

























































































hs P, , 2 | P, 
0.0 0.0 0.0 0.0 0.0 
0.2 2.19 2.20 2.04 214 
2.0 3.17 3.10 V1) 414 
40 1.90 3.90 a 86 390 
$4 4.55 4.55 4.$3 4.58 
S84 4.76 4.69 4.71 72 
- oe - ~-o- — a 
6.013 4.76 4.73 4.7] 476 
+ 7% 
6.015 4.88 4.82 479 4.90 
6.022 | 4.92 4.92 4.88 4.90 
6.040 6.04 § 92 5.90 § 90 
6.050 6.83 6.7) 67 673 
6 060 | 7.40 7.17 7.18 
6.070 7.50 7.33 7.46 
+ ——-+— ——--— 
6.080 7.3 7.40 7.46 
6.093 | 7.53 7.40 7% 











Note: P, through P, are designations of the sensors at different points in the combustion chamber 





The most interesting case was that of subjecting the 
engine to the very same pressure pulse, but in the middle 
of its operation, at which time there was still enough 
fuel in the combustion chamber to partially deform the 
nonstationary load and not permit catastrophic damage 
to une body. At the moment when the load was applied, 
the pressure in the combusuon chamber equaled 4.8 
MPa. It then increased smoothly to 7.5 MPa and 
subsequently did not significantly increase any further 
for 0.03 seconds until the body fractured (Table 2). The 
total time for which the solid-fuel-fired rocket engine 
operated after application of the load amounted to r 
approximately equal to 0.093 seconds. A significant 





increase in pressure (by 60 percent) (Figure 1) was 
evidently able to occur either because of the increase in 
combustion surface when fragments of the charge were 
chipped off or else because of disruption of the norm 
efflux of combustion products as a result of “blockage” 
of the critical cross section of the nozzle block by the 
chipped-off pieces of fuel. A similar situation occurred 
during tests of the SL-3 solid-fuel-fired rocket engine,'’ 
when several pieces of charge were torn off and carried 
away by the gas jet through the nozzle, which caused a 
sharp increase in pressure and fracture of the engine. 
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n 6} Mee 
Figure 1. Graph of P as a function of ¢: 1, nominal 
relationship; 2, relationship obtained in an 
experiment. (Y-axis represents P x 10°, Pa; X-axis 
represents time in seconds.) 


To estimate the time required for “blockage” of the 
critical cross section after applicauon of the load, 
we will proceed from the infinitesimality of the tme 
required for formation of the spalls, the sphericity of the 
chipped-off pieces (with a diameter commensurate with 
the size of the critical cross section R), and the equality 
of the velocity of the gas flow along the channel’s length 
and its average value V. From the equation for the fuel’s 
motion along the channel axis (where V is the velocity 
of a chipped-off piece), we obtain the time required for 
‘blockage,’ which has been designated 1, 


Pe LR 
w= 4( ice) ° 


where p, is the density of the fuel, p, is the density of the 
gas, L is the distance from the central part of the engine 
to the critical cross secuon, and C, 1s the aerodynamic 
drag given subsonic circulation. 


Calculations indicate that £, < 1°; consequently, blockage 
of the nozzle by the chipped-off pieces of fuel may 
occur. An increase in combustion surface when the 
charge is damaged is an equally likely reason for 
the increase in pressure inside the chamber, however. 
It follows from the relationship between the area of 
combustion S(t) and pressure inside the combustion 
chamber’, i.e., p/p(O) = S/S(0))S* (where v = 0.266 1s 
the index tn the pressure dependence of the combustion 
rate) that an increase in pressure to 7.68 MPe occurs on 
account of an increase in the area S$ in the amount of 
AS = 1,800 cm? (Table 3 presents data regarding the 
area of combustion given normal operation of a solid- 
fuel-fired rocket engine). A similar value is obtained 
for the increase in area if in view of the characteristic 
features of spalling in cylindrical charges,'”'* a crescent- 
shaped form of damage along the entire length of the 
combustion chamber is assumed and AS is esumated in 
accordance with the formula AS approximately equal 
to *R(2L) approximately equal to 1,880 cm*. 





Table 3. Dependence of the Area of Combustion on Thickness of the Fuel Arc 
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lcm S, cm2 lcm | S, cm? 
0 2,360 8.6 4,825 
1 2,776 9 4,817 
2 3,180 — 10 4,500 
3 3,558 11 —_ 4,216 
4 3,840 12 4.080 
§ 4,157 13 3,923 
6 4,000 14 4,935 
7 4,626 18 3,$27 
. | 4,810 16.2 2.325 
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The esumates made thus demonstrate the possibility 
f destrucuon of an engine on account of chipped-off 
pieces of fuel resulung from a nonstationary lateral load 
propagated in the charge by the shock. 


In conclusion, the authors express their gratitude to 
Z.G. Alpaidza, Kh.Kh. Ganeyev, V.I. Pakhunov, S.S. 
Solbodchikov, V.A. Tsvetkov, and A.A. Cheprunov for 
their great assistance in conducung the tests and for their 
valuable comments. 


References 


|. A.M. Sinyukov, L.I. Volkov, A.I. Lvov, A.M. 
Shishkevich, Ballisucheskaya raketa na tverdom toplive 
Solid-Puel-Fired Ballistic Missile), Voyenizdat MO 
SSSR, Moscow, 1972. 


2. V.1. Feodosyev, Osnovy tekhniki raketnogo po- 
ieta (Principels of Rocket Flight Engineering), Nauka, 
Moscow, 1979. 


4, A.V Karmuishin, E.D. Skurlatov, V.G. Startsev, V.A. 
Faldshteyn, Nestatsionarnaya aerouprugost tonkosten- 
nykh konstruktsty (Nonstauonary Aeroelastcity of 
Thin-Walled Structures), Mashinostroyeniye, Moscow, 
1982, 240 pages 


4. Ye.P. Velikhov, R.Z. Sagdeyev, A.A. Kokoshin, eds., 
Kosmicheskoye oruzhiye: dilemma bezopasnost (Space 
Weaponry: Safety Dilemma), Mir, Moscow, 1986, p 
| 82 


5S. S.P. Anisimov, A.M. Prokhorov, V.Ye. Fortov, 
USPEKHI FIZ. NAUK., Vol 142, No 3, 1984, p 395. 


6. Li. Sedov, Metody podobiya i razmernosti v 
mekhanike (Similarity and Dimensionality Methods in 
Mechanics), Nauka, Moscow, 1981, p 448. 


7. J.D. Achenbach, AIAA J., Vol 3, NO 6, 1965. 


§ W.E. Alzheimer, MJ. Forrestal, W.B. Murfin, ALAA 
J.. Vol 6, No 10, 1968. 


9 MI. Forrestal, W.E. Alzheimer, D.B. Longcope, 
AIAA J., Vol 8, No 9, 1970. 


10. M.I. Ugamov, V.A. Ivanov, B.V. Gulin, Prochnost, 
ustoychivost 1 dinamika obolochek s uprugim zapol- 
nitelem (Strength, Stability, and Dynamics of Shells 
With an Elastic Filler), Nauka, Moscow, 1977, p 331. 


11. MLL. Dgamov, V.A. Ivanov, B.V. Gulin, Raschet 
obolochek s uprugim zapolnitelem (Designing Shells 
With an Elastic Filler), Nauka, Moscow, 1987, p 260. 


12. VI. Pozhuyev, IZV. SSR. MMT, No 67, 1978, p 
106 


FBIS-UST-95-052 
19 December 1995 


13. A.V. Ostnk, I.B. Petrov, V.P. Petrovskiy, MATEM- 
ATICHESKOYE MODELIROVANIYE, Voi 2, No 8, 
1990, p 51. 


14. Dl. Bynum, L.U. Rastrelli, R.C. De Hart, AIAA J., 
Vol 2, No 2, 1964. 


1§. N.L. Coleburn, AIAA J., Vol 4, No 3, 1966. 


16. NJ. Murri, D.R. Curran, L. Seaman, Shock Waves 
in Condensed Matter, New York, AIP, 1982, p 460 


17. A.V. Oktrik, V.P. Petrovsky, Eighth Intern. Cont 
Equations of State, Elbrus, 1992, Abstracts, p 32. 


18. A.V. Oktrik, V.P. Petrovskiy, ZAPMTF, No 1, 1993, 
p 135. 19. MISSILE/SPACE DAILY, Vol 26, No 9, 
1967. 


Two-Dimensional Spalling of Cylindrical Polymer 
Structures 


957A1194A Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 4-10 


[Article by A.V. Ostrik and V.P. Petrovskiy, Central 
Physicotechical Institute, Russian Federation Ministry 
of Defense; manuscript received 22 Sep 94; UDC 
536.424:669.018.824) 


(FBIS Abstract] A method of calculating acoustic waves 
in cases where a filled infinitely long cylindrica! shell 
with a complexly shaped opening is under an axisym- 
metric load was proposed in a 1990 publication. The 
proposed method has now been generalized for the case 
where ignoring shearing waves in favor of the compres- 
sion waves would not be unwarranted. The mathemau- 
cai analysis focuses on the effect of a nonstationary load 
on a polymer-filled cylindrical shell that has an inner 
channel with a virtually random but rather smooth shape 
specified by a law governing the change in thickness of 
its filler at the angle coordinate 4(@). It is assumed that 
the conditions of plane deformation are satisfied and that 
Hook's law may validly be used to describe the poly- 
mer filler’s deformation properties. The equation system 
constructed to describe the elastic waves in the filler in 
a cylindrical coordinate system was derived on the ba- 
sis of equations proposed by Sneddon and Radek. The 
three-equation system derived to find the funcuons P 
(pressure), 4 o(the difference between the stress tensor 
components o, and ¥,), and 7, (the third component of 
the stress tensor) is a set of wave equations that can all 
be solved effectively by numerical methods. The equa- 
tion system's boundary conditions have been formulated 
under the assumption: ihat the shell holding the filler is 
elastic, the radial displacements at the shell-filler _nter- 
face are equivalent, and the inner boundary is free from 
stress. A numerical solution is obtained by representing 
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the sought functions in terms of Fourier series. The re- 
sultant linear equation system contains the amplitudes 
of the harmonics of the three unknown funcuions and 
is solved by the matrix sweep method. The results of 
the calculations of the waves’ propagation in the filler 
and their subsequent reflection from the channel surface 
were compared with data from experiments in which 
60-mm-long cylinders (outer diameter, 40 mm) made 
of a highly filled elastomer (with a measured density of 
1.8 g/cm” were subjected to a shock wave created by 
electrical explosion of type AS aluminum foil that was 
20 yam thick, 20 jam wide, and had a surface area of 
15 cm. In the experiments, a maximum voltage of 35 
kV was applied to the foil, which was placed in direct 
contact with the loading surface. Cylinders with cylin- 
dric and star-shaped inner channels were studied. The 
spalling of those cylinders with cylindrical channels had 
a characteristic crescent shape, and the amount of mate- 
rial actually damaged decreased as the diameter of the 
inner channel was decreased. An enurely different pat- 
tern of spalling was observed in cylinders containing 
star-shaped channels: The load angles recorded in the 
experiments with those cylinders equaled 30, 60, 90, and 
120 degrees, respectively. Areas of fracture connected 
with the development of cracks between the rays were 
observed along with the characteristic zones of spalling 
near the free surface. The presence of rays and their 
orientation with respect to the center of the load sig- 
nificantly affected the nature of the fractures observed. 
Unlike in the case of a cylindrical channel, the main 
portion of the rarefaction waves’ energy was concen- 
trated in the area between the star's rays, even when 
the said area was not directly under the load center. 
The concentrauon of energy in the space between the 
rays was attributed to interference of the waves reflected 
from the side walls of the rays and the inner portion o: 
the channel. The pressure isolines calculated by using 
the proposed mathematical model and the fracture zones 
predicted by the model as a function of the loading an- 
gle and the geometry of the inner channel were found 
to be in rather good agreement with those occurring 
in the actual experiments. It was therefore concluded 
that the proposed mathematical model may be used for 
qualitative, and in the case of some parameters, quan- 
titative description of the characteristic features of the 
wave process in cylindrical polymer structures with an 
inner channel having a complex geometry. Figures 3; 
references 14: 9 Russian, 5 Western. 


Conditions of Explosive Transformation of 
Explosives in Engineering Systems During the 
Modeling of Accidents 

957A1194B Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 18-23 


[Arucle by G.S. Smirnov, A.M. Khelemendin. \ M 
Gerasimov, A.K. Zhiyembetov, MD. Polyakoy, VA 
Rodionov, N.D. Manachinskiy, and I.P. Khabarov, A! 
Russian Scientific Research Institute of Experimenta 
Physics, Arzamas-16; manuscript received & Aug “4 
UDC 532.593] 


[FBIS Abstract} The response of explosive Charges | 
bullets/fragmentauon elements, shock waves, and heat 
pulses was studied in a series of tests involving ex- 
plosives based on tnnitrotoluene [TNT], TATB 
togen, and hexogen. The explosion parameters were 
recorded by filming, electromechanical, and x-ray meth- 
ods and by measuring its pressure, temperature, and air 
shock wave. In the experiments, an explosive desig- 
nated the acceptor explosive was subjected to the force 
of a shock wave generated by a shock wave generator 
or by an explosive charge designated the donor explo- 
sive. The range of initiating pressures from 1.0 to 4.0 
GPa, which is the characterisuc range for the typical 
anomalous effects to which explosives in engineering 
systems are subjected, was examined. The effects otf 
botn one- and two-dimensional loads on the acceptor 
explusives were studied. The experiments examining ex 
plosives’ response to anomalously high heat were con- 
ducted with explosive charges enclosed in sealed cham- 
bers. The combined effects of impact and fire were stud 
1ed in experiments in which plane and cylindrical gaps 
were create! in the explosive charge so as to simulate 
fracture of the explosive and the initiauon and devel- 
opment of convective combustion. The explosives were 
gnited by applying an electrical current to a nichrome 
spiral or by heaung the container holding the explosive 
in the flame of burning liquid fuel. According to the ex- 
periments, in most accident scenarios the anomalous!) 
high pressures and/or temperatures to which explosives 
are exposed are insufficient to excite a mode of nor- 
mal detonation in them but instead excite a nonstauon- 
ary low-speed explosive transformation. Purther devel- 
opment of the said process depends on the properties ot 
the given explosive, the dimensions of the charge, and 
the design features of the engineering system of which 
the explosive is a part. Depending on these factors, one 
of three outcomes proved to be possible: propagation 
Of a quasi-stationary low-speed explosive transforma 
tion throughout the entire charge, transformation of the 
low-speed explosive transformation process into a high- 
speed detonation process, and transformation of the low- 
speed transformation process into an attenuating explo- 
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sive transformation process in which part of the explo- 
sive is either scattered or burned up. In the latter case, 
the decomposition of the explosive in a specific stage of 
the process is compensated for by a pressure drop result- 
ing from efflux of the pase: into the resultant holes and 
expansion of the casing holding the explosive charge. 
In the case of one-dimensional shock wave loading, the 
region of propagation of the quasi-stationary low-speed 
explosive transformation of open explosive charges (i.e., 
charges not enclosed in a casing) is bounded on the bot- 
tom by the value of the critical diameter and on the top 
by the threshold diameter at which the transition from a 
low-speed explosive transformation to high-speed deto- 
nation occurs. Placing explosives in a metal casing may 
Stabilize the propagation of low-speed explosive trans- 
formation. The process of low-speed explosive trans- 
formauon culminating in high-speed detonation is the 
result of a shock wave. The quasi-stationary wave of 
low-speed explosive transformation is a complex con- 
sisting of a compression wave (precursor) and plastic 
wave corresponding to the process of explosive trans- 
formation. As a quasi-stationary wave of low-speed ex- 
plosive transformation moves throughout an explosive 
charge, its profile undergoes a process of amplitude-time 
evolution during which the compression wave and re- 
gion of chemical transformation converge as a result of 
the increase in intensity of energy release in the region 
of chemical transformation. When an acceptor explosive 
is subjected to two-dimensional loading, the transition 
from low-speed explosive transformation to high-speed 
detonation may occur at a point lying on a continuation 
of the axis of the effect in the depths of the explosive 
charge because of lateral relief. A low-speed explosive 
transformation process can progress into a process of 
high-speed detonation when waves are reflected fr ~ 
the metal shells surrounding explosives. These and ou 

findings of the experiments performed must be taken 
into amount when predicung the consequences of pos 
sible accidents triggering processes of low-speed explo- 
si\¢ transformation. Figures 6; references 6 (Western). 


Phenomenon of Dispersion of Powder Explosives 
Upon Pulse X-Irradiation 
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[Article by B.R. Gafarov, Ye.S. Yushkov, V.P. Yefre- 
mov, and Ye.B. Druzhkov, Chemical Technology Sci- 
entific Research Institute, Dzerzhinskiy; manuscript re- 
ceived 21 Sep 94; UDC 623.454.86] 


[FBIS Abstract] The phenomenon of dispersion of 
powdered octogen and ammonium perchlorate when 
they are subjected to a pulse of x-radiation lasting 
approximately a nanosecond was studied. The octogen 
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crystals ranged in size from 50 to 1,S00 jm, and 
the ammonium perchlorate granules ranged in size 
from 30 to 1.000 jm. Both finely dispersed (150-300 
ym) and coarsely dispersed (500-700 jzm) fractions 
of each of the two explosives were studied. Their x- 
radiation pulses lasted only a few nanoseconds each, 
and their energy contribution ranged from 20 to 40 J/ 
g. During the experiments, 4-mm layers (SO-g portions) 
of explosive were poured into acrylic plastic containers 
and covered with a 3-mm-thick sheet of acrylic plastic, 
which filtered out the soft component of the x-radiation 
and equalized the energy release throughout the layer. 
The container was irradiated in the open air. Individual 
granules of explosive and pressed polycrystalline wafers 
of potassium bromide and cerium nitrate (thickness, 
0.25-0.5 mm) were also irradiated under the same 
conditions. Irradiation with a pulse electron accelerator 
and static irradiation on a gamma unit (“Co) were 
also performed. The irradiated and control specimens 
were subjected to a comprehensive physicochemical 
analysis that included optical and scanning electron 
microscope studies, differential thermai analysis, x-ray 
crystallographic analysis, and x-ray phase analysis. The 
most pronounced effect observed, 1.e., dispersion of the 
large fracuons of the explosive powder, was induced 
by the x-ray pulse irradiation. A particle size fraction 
threshold of the onset of fracture, 1.e., 200-300 jim, 
was discovered. In addition, irradiation reduced the 
particles’ strength. The octogen crystals became cloudy; 
the radiation-chemical damage observed after they had 
been subjected to pulse x-irradiation was also evident 
after they had been irradiated under the other conditions 
Studied. The alterations in their remaining parameters 
remained within the bounds of the errors associated 
with the various measurement methods. All of the 
pressed wafers of salts were totally broken down into 
a powder with a particle size below 5 ym. The largest 
individual granules were also broken down into very 
fine fragments. Three mechanisms of brittle fracture that 
may possibly occur during so highly intense an effect as 
pulse x-irradiation were discussed: 1) thermomechanical 
Stresses aS a result of uneven heating; 2) rapid gas 
evolution as a result of radiolysis; and 3) shock wave- 
induced brittle fracture. Table 1; references 3: 2 Russian, 
1 Western. 


Explosive Tubular Accelerator for High-Speeding 
Throwing 


957A1194D Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 27-32 


[Article by A.Yu. Dolgoborodov, Chemical Physics 
Institute imeni N.N. Semenov, Russian Academy of 
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Sciences, Moscow; manuscript received 17 Aug 94; 
UDC 534.222.2] 


[FBIS Abstract] An explosive tubular accelerator has 
been proposed for use in processes of throwing metal 
disks at high speeds. The proposed accelerator makes 
it possible to increase throwing speeds and reduce the 
weight of the explosive charge in relation to the weight 
of the body being thrown. The accelerator includes a 
hollow cylindrical high-power explosive charge whose 
core is filled with a composite charge consisting of 
a lightweight explosive with a high hydrogen content 
and a heavy, dense explosive located directly after the 
initiating system in front of the lightweight explosive. 
The metal disk that is to be thrown is accelerated in 
a special barrel. In existing explosive tubular accelera- 
tors the pulse is communicated to the thrown body by 
compressed gas formed from the foamed polymer ma- 
terial that fills the accelerator’s core. In the proposed 
explosive tubular accelerator, the thrown body is ac- 
celerated by a fundamentaliy different method. Specif- 
ically, it is accelerated by filling the accelerator’s core 
with an explosive charge with a relatively high hydro- 
gen content and thereby increasing the throwing speed 
while reducing the amount of explosive used. On the 
basis of previously reported studies examining the opti- 
mum conditions for forming a Mach detonation wave in 
explosive compounds inside hollow cylinders of more 
powerful explosives, a 50:50 mixture of octogen and 
tungsten was placed in the front part of the core for 
additional supercompression of the products generated 
when the nitromethane explodes. The products of ni- 
tromethane’s supercompressed detonation serve to ac- 
celerate the metal disk being thrown. In tests, the force 
generated by the mixture of lightweight and heavy ex- 
plosives and the dense octogen and tungsten charge 
proved sufficient to accelerate metal disks with a thick- 
ness of 2 mm and diameter of 8 mm. The duralumin 
disks accelerated by the explosive tubular accelerator 
reached a velocity of 8 km/s, and the steel disks reached 
7 km/s. Additional research is planned to further in- 
crease the new accelerator’s throwing speed by using 
lightweight explosives containing even higher levels of 
hydrogen and by using more powerful compounds for 
the heavy and high-power explosive charges. Figures 2; 
references 9: 7 Russian, 2 Western. 


Temperatures of Attainable Superheating of 
Energy-Storing Materials Under Conditions of 
Thermal Decomposition 
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[Article by O.F. Shienskiy and V.M. Aristov, Rus- 
sian Chemical Technology University imeni D.I. 


Mendeleyev, Moscow; manuscript received 12 Apr 94 
UDC 678.664] 


(FBIS Abstract] The method of contact healing was used 
to study the kinetics of the thermal decomposition of se 
lected condensed explosives. Study specimens of tused 
explosives ‘vere applied to metal substrates as thin films 
The thickness of the films, 1.e., 2 to 7 jam, made 11 pos 
sible to quickly heat the study specimens to the temper- 
ature of the metal substrate. Study specimens in powder 
form with particle sizes less than 0.01 mm were placed 
in a melt of tin or in PES-S liquid silicone. The speci- 
mens were placed on a metal substrate with a caulked 
thermocouple. Measurement of the time required for 
completion o: the thermal decomposition reaction [1 | 
was based on cessation of the release of bubbles ot 
the gasiform reaction products. The experiment results 
were represented in the form of graphs of lgi l/r.) = 
fU\/T). The experiments revealed a previously unknown 
property of condensed explosives, 1.e.. the fact that the 
temperature range of the thermal decompositon process 
is bounded by the limiting temperatures. The activation 
energy values obtained for the explosives trinitrotoluene 
[TNT], ammonium nitrate, sodium azide, TATB, nitr 
cellulose, ectogen, ammonium bichromate, lead azide 
and petroleum from the Alanin oil field corresponded 
to published activation energies, and the graphs plotied 
for Ig(1/t.,) = f(1/T) were in satisfactory agreement with 
the empirical equation previously proposed for linear 
polymers. in the case of TNT, the limiting temperature 
of the thermal decomposition process [7)] was close | 
its boiling point (320-340°C). The value of 7. found for 
ammonium hichromate corresponded to the temperature 
in its thermal Cecomposition reaction front upon excita- 
tion of low-speed detonation (300°C). The iemperature 
reached by nitrocellulose during thermal decomposition 
corresponded to that published elsewhere. It was con 
cluded that in a physical sense, the temperature 7, char- 
acterizes the intensification of homogeneous nucleation 
center formation in individual atoms of the crystalline 
lattice and in vacancies. The sharp bend observed in 
the graph of lg(1/t.,) = fl1/T) was explained by the fact 
that the rate of exponential increase in the frequency of 
formation of homogeneous nucleation centers and the 
number of heterogeneous micleation centers in disloca- 
tions and vacancies is steeper than the increase in the 
Arrhenius exponent of the chemical reaction. Possible 
reasons for the observed violation of the Arrhenius rela- 
tionship were offered. Table 1; references 10: 7 Russian, 
3 Western. 
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Zherebtsov, N.I. Leonova, Z.V. Surayeva, and A.V. 
Timonkina, All-Russian Scientific Research Institute of 
Experimental Physics, Arzamas-1!6; manuscript received 
15 Aug 94; UDC 536.71} 


[FBIS Abstract] Sets of problem-oriented applications 
packages have been developed at the Mathematics De- 
partment of the All-Russian Scientific Research Institute 
of Experimental Physics (VNIIEF] for solving various 
scientific-technical and engineering problems. Each of 
the programs permits the solution of a specialized equa- 
tion system. The equation systems used in the programs 
are generally supplemented (closed) with equations de- 
scribing the thermophysical properties of materials used 
in structures being designed. Among the available equa- 
ion sets are characteristic equations linking three ther- 
modynamic quantities (most often density, pressure, and 
temperature or density, pressure. and energy) designed 
for use in cases of high pressures and temperatures. Ini- 
tially, each specific set of programs had its own library 
of characteristic equations. Eventually, several versions 
of a characteristic equation library were developed for 
maintaining all of the individual program packages de- 
signed for solving one-dimensional problems. The latest 
version, which is known as the URS-YeS-OK (char- 
acterisuc equation library for solving one-dimensional 
problems on YeS computers), currently contains ap- 
proximately 20 applications programs combined into 
a general service (utility) system and makes it pos- 
sible to fill requests for calculating approximately 90 
different quantities. More recently, a team of mathe- 
maticians from the All-Russian Scientific Research In- 
stitute of Experimental Physics and All-Russian Scien- 
tific Research Institute of Engineering Physics has de- 
veloped an improved equation library called the Gen- 
eral Pund of Characteristic Equations (URS-OF] that 
contains the equations required to solve virtually all 
applied problems written in algorithmic language. To- 
gether, the URS-YeS-OK and URS-OF contain more 
than |,000 characteristic equations for various materials 
and make it possible to solve a wide range of classes 
of problems, including problems dealing with the fol- 
lowing: local, wide-range, and global states; problems 
involving tabular, tabular-analyuc, and analytic data; 
and problems involving pure elements, alloys, chem- 
ical compounds. mixtures, minerals, condensed mate- 
nals, and gases. Work in the following areas is now 
under way at the Mathematics Department of the VNI- 
IEF: maintenance of the characteristic equations cur- 
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renuy existing in the URS-YeS-OK and URS-OF; mod- 
ification of existing characterisuc equauions (improving 
their precision, taking new effects into consideration, 
expanding the range of describable states, etc.); and 
analysis of models describing the properties of mate- 
rials and development of new characteristic equations. 
A program package called the Characterisuc Equation 
Center [TSURS] that is designed for use by interested 
specialists has been created for studying characteristic 
equauions and consists of three sections as follows: li- 
brary of programs for solving various standard prob- 
lems; library of analytic characteristic equations and 
sets of constants for them; and compendium of vari- 
ous experimental and theoretical data regarding materi- 
als’ thermophysical properties. The components of the 
TsURS have been briefly described along with two ad- 
ditional program packages that are being developed by 
the VNITEF for approximating tabular data with splines. 
References 26 (Russian). 


Calculating Thermodynamic Functions in a Liquid- 
Vapor Two-Phase Region 
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(FBIS Abstract] A method is proposed for calculating 
thermodynamic functions in a two-phase liquid-vapor 
region by using analytic relauonships to specify the in- 
dependent variables density and specific energy. The 
use of analytic relationships rather than iteration algo- 
rithms to specify density and specific energy 1s said to 
Significantly reduce the amount of computer time re- 
quired to calculaie the thermodynamic funcuons. The 
proposed calculation method is not ted to any spe- 
cific type of characteristic equation. All that is required 
to calculate thermodynamic functions and their deriva- 
lives is information regarding the state of the material 
at the boundary of the two-phase region. The required 
parameters at the said region are assumed to be known 
and representable in the form of polynomial (cubic and 
quadratic) splines. First, the states of copper in the two- 
phase liquid-vapor region were calculated on the ba- 
sis of tables constructed from splines obtained by us- 
ing the proposed calculation method. The said states 
were then recalculated by using the conventional itera- 
tion algorithm for calculating thermodynamic funcuons 
in a two-phase region. When a YeS computer was used 
for both sets of calculations, the proposed method re 
quired only about half the time required with the itera- 
tion method. It was noted that the proposed calculation 
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method only requires the use of quadratic splines when 
calculating certain specific functions. In most cases, the 
thermodynamic funcuons describing states at a liquid- 
vapor interface may be represented by any relationships 
whose derivatives can be calculated, including funcuons 
derived from physical models. Figures 2; references § 
(Russian). 


Model of Wide-Range Characteristic Equations of 
Polymeric Materials With High Energy Densities 


957A1194H Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 47-52 


[Article by I.V. Lomonosov, V.Ye. Fortov, and K.V. 
Khishchenko, Thermal Physics of Pulse Processes Sci- 
entific Research Center, Russian Academy of Sciences, 
Moscow; manuscript received 6 Oct 94; UDC 536.2] 


(FBIS Abstract} A caloric model of characteristic equa- 
tions has been developed that makes it possible to ef- 
fectively describe the properties of polymer materials 
within a wide range of densities and pressures, includ- 
ing in condensed and quasi-gaseous phases. The pro- 
posed model is a development of a previously proposed 
approach to constructing characteristic equations of ma- 
terials that are complex from the standpoints of chemi- 
cal composition and physical properties. The proposed 
model is specified in the following generalized form: 
PV, E) = pV) + [(T(V, EVVUE -E(V))), where E(V) 
and p.(V) = -dE /dV are the elastic components of the en- 
ergy and pressure when 7 = 0 K, and the coefficient I'(V, 
E) specifies the contribution of thermal components to 
the characteristic equation. The proposed model is used 
to calculate the thermodynamic characteristics of poly- 
ethylene, epoxy resin, and polyamide. The results are 
compared with published results of experiments involv- 
ing the said polymeric materials and high energy densi- 
ties. The two sets of results were found to be consistent 
with one another throughout the entire phase diagram of 
each of the materials studied. The proposed model was 
recommended for use in numerical modeling of nonsta- 
tionary hydrodynamic processes. Figures 3; references 
10: 4 Russian, 6 Western. 


Selected Features of the Effect of Ultradeep 
Penetration of Powder Particles Into a Metal 
Eaclosure 


957A11941 Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 53-58 


[Article by S.K. Andiievko, §.M. Usherenko, and V.A. 
Shilkin, Pulse Processes Scienufic Research Institute, 
Byelarus Republic Powder Metallurgy Scientific Pro- 


duction Associauon, Minsk; manuscript received 8 Aug 
94; UDC 532.593) 


[FBIS Abstract) A series of experiments were conducted 
to examine the processes occurring upon the interaction 
of dense (up to 5 g/cm’) high-speed (up to 3 km/s) 
flows of powder particles within a metallic enclosure. 
The experiments revealed an effect that is fundamentally 
new for the given region of interacuon. Specifically, it 
was discovered that a noticeable percentage of parti- 
cles of the flow (generally from 0.1 to | percent) pene- 
trate the metallic enclosures to depths as much as | ,000 
Or more times their initial diameter (caliber), whereas 
under normal conditions, the potenual energy of any 
particle of the flow is sufficient to cause the said parti- 
Cles to penetrate the enclosure to a depth not exceeding 
50 calibers. The newly observed effect was termed ul- 
tradeep penetration and was studied theoretically and 
experimentally. The experiments established that ultra- 
deep penetration occurs only under conditions of in- 
teraction of the metal enclosures and the dense flows 
of powder particles. Transmission electron microscopy 
and x-ray microanalysis demonstrated that penetration 
of the particles into the enclosure material is accom- 
panied by well-developed intensive plastic deformation 
with strong heating of the glide planes directly adjacent 
to the penetrating particle. Numerous studies revealed 
that the channels formed in the enclosure by the part- 
cles as they penetrate it collapse after they have passed 
through. As the particle penetrate the metal enclosure, 
their mass is significantly reduced on account of entrain- 
ment of a portion of the mass of their boundary layer 
as they interact with the enclosure. The trajectories of 
the particles’ motion in the material of the enclosure 
are not always a Straight line. Experiments in which Ti, 
Fe, Cu, Al, Pb, and various types of steel functioned as 
the target and nitrides, borides, silicides, and powders 
of such metals as tungsten, iron, and lead functioned as 
the penetrating powder indicated that ultradeep penetra- 
tion largely depends on the chemical compositions of 
the interacting objects. Figures 5; references 2 : 1 Rus- 
sian, | Western. 


Compaction of Ultradisperse Diamonds in a Weak 
Shock Wave 


957A1194J Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 59-61 


[Article by E.E. Lin, V.A. Medvedkin, and S.A 
Novikov, All-Russian Scientific Research Institute 
of Experimental Physics, Arzamas-16; manuscript 
received 8 Aug 94; UDC 532.593} 


[FBIS Abstract] The compaction of ultradisperse dia- 
monds in a weak shock wave was examined in a series 
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of experiments in which steel pistons with the appro- 
priate linear dimensions were used to accelerate vials 
filled with ultradisperse diamond powder to velocities 
of 500 +/-50 m/s. Analytical esumates confirmed that a 
pressure pulse with an amplitude in the range of weak 
shock waves had indeed been created in the system of 
ultradisperse diamond clusters studied. Analyses of the 
study samples based on analyucal chemistry techniques, 
electron and optical microscopy, and x-ray diagnosis 
established that after the shock wave loading, the pow- 
der consists primarily of polycrystalline diamond par- 
ucles with a cubic modification of their crystal lattice. 
The region of coherent scattering of x-radiation on the 
resultant parucles corresponded to the average size of 
the ultradisperse diamond. The histogram of the parti- 
cles of diamond powder after dynamic compaction of 
the ultradisperse diamond could be approximated by a 
lognormal distribution. Individual, comparatively large 
grains up to 600 ym in diameter were detected in the 
study specimens. The experiments thus confirmed that 
when ultradisperse diamond powder is subjected to the 
impact of a weak shock wave, the average size of the 
diamond particles increases by more than two orders of 
magnitude. It hypothesized that the said effect is in all 
likelihood due to background excitations of the crystal 
lattice of the clusters of ultradisperse diamond in the 
weak shock wave. Figures 2; references 7 (Russian). 


Macrokinetic Features of the Occurrence of a 
Two-Way 6-a Phase Transformation in 
46-Plutonium With Threshold Load Pulses 


957A1194K Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan95 Vol 14 No 1, pp 62-68 


(Article by Ye.A. Kozlov, L.F. Timofeyeva, and A.K. 
Muzyrya, All-Russian Scientific Research Institute of 
Thermai Physics, Snezhinsk (Chelyabinsk-70), and All- 
Russian Scientific Research Institute of Inorganic Mate- 
nals imeni Academician A.A. Bochvar, Moscow; man- 
uscript received 22 Sep 94; UDC 532.593} 


[FBIS Abstract] Static experiments have established that 
the pressure of the initiation of a 4 — a phase tran- 
sition in 4-plutonium depends on the concentration of 
alloy-forming elements (aluminum, gallium, etc.) in the 
4-plutonium alloy, whereas the presence of shearing 
Stresses facilitates the occurrence of the said transfor- 
mation in accordance with a martensite mechanism. In 
experiments wherein 4-plutonium has been subjected to 
the effects of weak shock waves, drastic variation in the 
velocity of the plastic wave has been observed. In view 
of these facts, a series of experiments were conducted 
to study the formation of a complex wave configuration 
in thin layers of a 4-phase alloy of plutonium and gal- 
lium when the alloy was subjective to explosive loading 
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in the vicinity of the two-way 4-a phase transforma- 
tion. The primary objective of the experiments was to 
identify macrokinetic laws governing the occurrence o! 
two-way 4-a phase transformauions. The experiments 
were performed on three samples of one and the same 
plutonium-gallium alloy in the form of a solid solu- 
tion based on 4-phase plutonium. The total amount of 
impurities in the alloy did not exceed 0.20 percent by 
weight. The samples were wedge shaped and cut from 
parallelepipeds measuring 10 x 30 x 70 mm. The an- 
gle at the vertex of each sample equaled 12°) 00’ +/ 
- 5’. The experiments were performed in a specially 
constructed sealed unit permitting continuous record- 
ing of the wave processes occurring in the study sam- 
ples throughout the pulse loading process. When the ¢- 
plutonium alloys were subjected to the effects of weak 
shock waves generated by detonation of an explosive 
charge, a three-wave configuration developed. Unlike 
the front of the first wave, the fronts of the second and 
third waves were drawn out and underwent a gradu?! 
increase in mass velocity. A time delay in formation of 
the main plastic wave in the material at threshold load 
pulses was established along with a rapid slowing of the 
increase in the loading pulse’s amplitude and duration 
Scanning electron microscopy studies of intact samples 
of the 4-plutonium alloy confirmed that the observed 
phenomenon is linked to nonuniform distribution of the 
alloy-forming element throughout the grain and to the 
occurrence of a 6 — a phase transformation at the lower 
stress levels [7,,] in those segments of the grains that 
are relatively poor in alloy-forming element. It was fur- 
ther discovered that as the loading pulse’s amplitude and 
duration increase to levels sufficient for the occurrence 
of a phase transformation not just in those segments of 
6-plutonium grains adjacent to the grain boundaries but 
throughout the bulk of the grains as well, the delay in 
formation of the main plastic wave decreases rapidly 
and drops to zero when o,, 2 o°,,, which corresponds 
to the transition from a three-wave (elastic precursor, 
fuzzy phase precursor, and main plastic wave) to a (wo- 
wave configuration in é-plutonium. Figures 4, table |; 
references 17: 6 Russian, 11 Western. 


Behavior of a-Quartz Under High Dynamic and 
Static Pressures: New Results and Ideas 


957A1194L Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 69-74 


[Article by Yu.N. Zhugin, All-Russian Scientific 
Research Institute of Technical Physics, Snezhinsk 
(Chelyabinsk-70); manuscript received 21 Sep 94; 
UDC 539.89:549.51] 


[FBIS Abstract] The notion of the absence of non- 
Stauonary splitting of the shock wave front in quartz 
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throughout the entire region of phase mixing beginning 
at 12-14 GPa and lasting an anomalously long time from 
the standpoint of the differences in the stress (i.e., pres- 
sure) level between the beginning and completion of the 
transformation of quartz into stishovite (more than 20 
GPa) had been established by 1990. Many researchers 
also believed that at stress levels (pressures) exceeding 
the dynamic elastic limit, the quartz behind the shock 
wave front did not possess any shearing strength. The 
year 1980 marked the beginning of a shift in think- 
ing regarding the effects of high levels of stress and 
pressure on quartz as fundamentally new research find- 
ings emerged from dynamic experiments, static exper- 
iments, and molecular dynamics calculations. Publica- 
tions summarizing those experiments and calculations 
were reviewed and analyzed, and their findings were 
compared. Specifically, when a hypothetical shock adia- 
bat was determined for quartz based on published results 
of laboratory experiments involving samples of quartz 
subjected to one-time shock wave loading, its position 
in pressure-volume [P-V] coordinates turned out to be 
fundamentally different from the accepted shock adiabat 
for quartz. The hypothetical shock adiabat virtually co- 
incided with the curve plotted for one-time shock wave 
compression of quartz as estimated in macrokinetics in- 
volving mixtures of quartz with paraffin and fluoroplas- 
tic. In the range of pressures from 23 to 35 GPa, both 
curves in P-V coordinates were located below the ray 
corresponding to the critical pressure of 23 GPa, which 
Causes splitting of the shock front in the said pressure 
range. When the results of shock wave and static exper- 
iments were compared, the similarities of the pressures 
corresponding to the following proved striking: the be- 
ginning of the gradual anomalous compression of quartz 
in dynamic experiments (8-14 GPa) and the beginning 
of its gradual progression to an amorphous state (12.5 
GPa) in static experiments; the beginning of its more 
abrupt compression (20-23 GPa) and the beginning of 
its more noticeable progression to an amorphous state 
(20-25 GPa); and the beginning of steep legs in the hy- 
pothetical shock adiabat and curve plotted for one-time 
shock compression of quartz (31 GPa) and the poit.t at 
which quartz becomes completely amorphous (30 GPa). 
When the dynamic P-V curve was compared with the 
results of molecular dynamics calculations, it was dis- 
covered that the experimentally determined critical pres- 
sure of 20-23 GPa virtually coincides with the pressure 
at which quartz’s volume abruptly decreases by approx- 
imately 15 percent. In addition, the pressure level of 
approximately 30 GPa that is characteristic for dynam- 
ics experiments was found to absent in the molecular 
dynamics calculations. It was therefore concluded that 
given the present state of research, there are sufficient 
grounds io view the position of the curve of one-time 
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shock compression and the hypothetical shock adiabat 
of quartz in P-V coordinates as just a second approxima- 
tion. Figures 2; references 35: 13 Russian, 23 Western. 


Effect of the Initial Unevenness of a Material’s 
Density on Shock Wave Motion 


957A1194M Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 75-88 


[Article by V.A. Shcherbakov, All-Russian Scientific 
Research Institute of Experimental Physics; manuscript 
received 1 Sep 94; UDC 532.593] 


(FBIS Abstract) The problem of the motion of a 
plane, initially stationary shock wave through a solid 
body having an uneven density was solved. The shock 
wave was assumed strong enough that the strength and 
viscosity of the solid body could be ignored yet weak 
enough that the effects of heat conduction could also be 
ignored. It was further assumed that viscosity and heat 
conduction could only affect the structure of the shock 
wave front, which was of no interest in the problem at 
hand. The analysis focused on the curving of the plane 
shock wave as it was propagated through the solid body. 
In an approximation of density perturbations an analytic 
solution was obtained that established the fact that given 
a shock wave with a specified intensity, the effect of 
the unevenness of the density of the material through 
which that shock wave is propagated on the shock 
wave’s curvature during propagation will be attenuated. 
The case of porous polystyrene (foamed plastic) was 
considered as a specific example. Figures 2; references 
5 (Russian). 


Numerical Modeling of the Dynamics of Shock 
Wave Propagation in Composite Materials 


957A1194N Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Voi 14 No 1, pp 89-96 


[Article by M.-F. Ivanov, A.N. Parshikov, and V.A. Gal- 
burt, Thermal Physics of Pulse Effects Scientific Re- 
search Center, Russian Academy of Science’ Moscow; 
manuscript received 6 Oct 94; UDC 532.529:539 89] 


[FBIS Abstract] The dynamics of the propagation of 
shock waves in composite materials was examined by 
way of numerical modeling of the propagation of a 
compression pulse across fibers of a unidirectional rein- 
forced composite metal (an aluminum matrix with tung- 
sten fibers). Propagation of a shock wave through the 
composite metal was modeled in terms of the planar 
two-dimensional flow of an elastoplastic medium whose 
dynamics were represented in Lagrangian coordinates. 
The composite’s microstructure was given explicit con- 
sideration, and its elastoplastic behavior was described 
in terms of equations for the components of the elastic 
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stress deviator. The condition of Mises fluidity was used 
to determine the moment of the material's transition 
from an elastic state to a plastic state. The equations con- 
structed for the elastoplastic medium’s dynamics were 
solved by the method of finite differences on a La- 
grangian triangular net. The dependence of the scatter of 
the pulse transmitted through the layer of composite on 
the coefficient of three-dimensional reinforcement and 
the diameter and material of the reinforcing fibers was 
considered. The geometric scatter of the stress (pres- 
sure) pulse was estimated through the time required for 
the pulse to be transmitted to the rigid wall. It was de- 
termined that the scatter of the pressure pulse increases 
as the duration of the incident pulse decreases (which 
is consistent with the conclusions of the linear theory), 
and the time required for the pulse to be transmitted 
to the rigid wall increases as the strength properties of 
the composite decrease. The calculated scatter increased 
from 1.67 to 2.16 when stratification of the fibers and 
matrix in the case of an absolutely © astic composite 
was taken into consideration. In the case of elastoplastic 
deformations, stratification was not found to have any 
effect on scatter of the pressure pulse. Pulses passing 
through the composite with a three-dimensional rein- 
forcement coefficient of 0.4 were found to experience 
the maximum amount of scatter. The said result was 
found to be consistent with data from experiments in 
which a unidirectional composite with an epoxy matrix 
and aluminum oxide fibers (30 zm in diameter) was sub- 
jected to shock wave loading. Next, a relaxation model 
of a homogenous medium possessing dispersing and dis- 
sipative properties of the study composite was formu- 
lated. The validity of the approximation used to describe 
the dissipative and dispersing properties of the unidi- 
recuonal composite subjected to intensive pulse load- 
ing was confirmed by comparing the results of calcula- 
tions performed by using the one-dimensional relaxation 
model with results obtained by using the aforementioned 
two-dimensional model. Figures 3; references 7: 5 Rus- 
sian, 2 Western. 


Generation of Strong Shock Waves by Soft X- 
Radiation of Z-Pinch Plasma 


957411940 Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 96-100 


[Article by Ye.V. Grabovskiy, O.Yu. Vorobyev, K.S. 
Dyabilin, M. Ye. Lebedev, V.P. Smirnov, and V.Ye. For- 
tov, Innovative and Thermonuclear Research Institute, 
Troitsk, and Thermal Physics of Pulse Effects Scientific 
Research Center, Moscow; manuscript received 8 Aug 
94; UDC 532.593] 


[FBIS Abstract] A pulse of soft x-radiation was used 
to generate shock waves with amplitude pressures up 
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to 3 Mbar in targets of condensed aluminum and lead. 
The x-radiation used in the experimen’ had a power 
level of up to 2 TW/cm? and was creaied hy dynamic 
compression and irreversible heating of plasma in two- 
Stage Z-pinch cascade targets on an ANGARA-S-] 
unit. Using the Z-pinch plasma-based method made 1t 
possible to create pulsed x-radiation with a duration 
and a strength an order of magnitude greaer than 
the duration and strength of pulses created by the 
units converting laser radiation into x-radiation that are 
generally used in similar experiments. To generate the 
radiating plasma, the researchers used an inner liner 
in the form of a hollow cylinder consisting of a low- 
density agar-agar framework into which molybdenum 
had been implanted (so that the density of the cylinder 
was less than 10 mg/cm’ and its total weight amounted 
to less than 200 yg). The liner had an outer diameter 
of 4 mm, wall thickness of 0.2 mm, and height of 10 
mm. An ultrasonic annular stream of xenon weighing 
150 yg served as the outer liner. A 3.5-MA pulse 
current was passed through the outer liner. The Z-pinch 
plasma emitted soft x-radiation with a temperature on 
the order of 60-120 eV. The radiation hit a plane target 
that consisted of layers of aluminum (16 jim thick) 
and lead (80-200 jm thick) and that was located | 
mm above the inner liner. The optical basis method 
was used to record the shock wave’'s velocity. Quartz- 
polymer fiber optic communication lines were used to 
lead the radiation out of the experimental assembly. 
The average rate of propagation of the shock wave 
equaled 7.3 +/- 0.6 km/s on an 80-ym base and 4.6 
+/- 0.3 km/s on a 200-ym base, which according to 
the shock adiabat of lead, corresponds to average shock 
compression pressures of 3 and 0.9 Mbar, respectively. 
The propagation rates are also determined by numerical 
modeling under the assumption that the plasma had a 
brightness temperature of 80 eV. The results of the 
experiments and numerical modeling were cc isistent 
within the bounds of the measurement error. Figures 
3; references 16: 10 Russian, 6 Western. 


Examples and Prospects of Using Spherically 
Convergent Shock Waves To Study Physicochemical 
Transformations of Minerals and Mineral Ores 


957A1194P Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 108-119 


[Article by Ye.A. Kozlov, Yu.N. Zhugin, B.V. Litvinov, 
V.N. Kozlovskiy, and G.V. Kovalenko, All-Russian 
Scientific Research Institute of Technical Physics, 
Snezhinsk (Chelyabinsk-70); manuscript received 22 
Sep 94; UDC 5494539. 89+5 36.4] 


[FBIS Abstract] Equatorial secuons of selected miner- 
als and mineral ores that had been subjected to loading 
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by a spherically convergent shock wave were exam- 
ined and compared with intact spherical specimens of 
the same minerals/mineral ores. The following mineral 
ores were studied: quartzite (from the Pervouralsk de- 
posit), granosyenite (from Kazakhstan), marble (from 
the Gazgan deposit), serpentinite (from the Urals), and 
forsterite (from the Kola Peninsula). The study samples 
all measured 48 mm in diameter at the start of the ex- 
periments. In a vacuum (10° torr) they were sealed into 
pressure casings made of 12Kh18Nil0T steel with a to- 
tal shell thickness of 8 mm and outer diameter of 64 
mm. A sphere and two control samples were cut from 
each lump. Each of the sealed balls was then loaded 
with spherically convergent shock waves generated by 
detonating a spherical layer of explosive consisting of 
trinitrotoluene and hexogen (3:7). The explosive was 
triggered at a radius of 40 mm. During and after the ex- 
plosions, the pressure casings maintained their strength 
and seal, thereby excluding the possibility of contami- 
nation of the study materials. All of the specimens sub- 
jected to the spherically convergent shock waves and 
the control specimens were subjected to gammascopy 
in two mutually perpendicular directions and were then 
cut along a plane in the direction perpendicular to con- 
tact with the explosive. The cut sections were subjected 
to visual observation and examination with an electron 
microscope. According to the estimates presented, the 
pressures created in the various study samples at the 
front of the convergent shock wave ranged from ap- 
proximately 20 GPa on the samples’ surfaces to ap- 
proximately 10° GPa and temperatures of approximately 
10° K in a l-mm radius given a loading pulse dura- 
tion of approximately 0.5 ys. Of special note was fact 
that the shock wave-loaded samples of marble, serpen- 
tinite, and forsterite contained a central vacuum cavity, 
whereas the shock wave-loaded samples of quartzite, 
granosyenite, and graphite did not. This difference was 
hypothesized to be related to a change in the nature 
of energy accumulation at the front of the spherically 
convergent shock wave in quartzite, granosyenite, and 
graphite as a result of their distinct polymorphous trans- 
formations in the shock wave (which were accompanied 
by big changes in their specific volumes) and to a pos- 
sible filling of small cavities with solidified melt under 
the conditions of the experiments performed. The load- 
ing pulses reached in the intact mineral samples, which 
were considered to be “geothermal barometers,” were 
said to be known for certain to be sufficient not only 
for Ux ~ccurrence of numerous polymorphous transfor- 
mations but also for the occurrence of melting, vapor- 
ization, and (in a number of cases) thermal decomposi- 
uon into simpler chemical compounds. The uniqueness 
of the information presented as compared with informa- 
tion presented by other researchers is said to lie in the 
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fact that it details the appearance of all possible stages 
of the impact metamorphism of mineral ores at vari- 
Ous points along their radii. Figures 6; references 22: 
13 Russian, 9 Western. 


Amplitude and Duration of Shor: Axisymmetric 
Convergent Mechanical Pulse 


957A1194Q Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 119-123 


[Article by B.R. Gafarov and V.P. Yefremov, Chemica! 
Technology Scientific Research Institute, Dzerzhinskiy; 
manuscript received 29 Sep 94; UDC 532.598:678.01 | 


{[FBIS Abstract] Formulas were derived for the ampli- 
tude and duration of a short axisymmetric triangular 
shock wave converging to the axis of a thick cylinder 
of a material with a linear shock adiabat. The problem 
of finding the amplitude and duration of such a wave 
was formulated in an axisymmetric cylindric geometry 
in a hydrodynamic approximation. The correctness of 
the proposed mathematical model was verified by per- 
forming a series of test calcu'ations and comparing them 
with the results of an experiment in which a Plexiglas 
cylinder with an outer diameter of 40 mm was subjected 
to the force generated by the electrical explosion of foil 
The comparison established that 1) taking the amplitude 
dependence of the wave velocity into account results 
in amplitude and duration values that differ severalfold 
from the values obtained with an acoustic approxima- 
tion and 2) focusing of the pulse will only be observed 
when the pulse’s length becomes comparable with the 
distance to the axis. It was concluded that taking non- 
linearity into account, even if only in a first approxi- 
matio”, is extremely importart when calculating short 
high-i ‘ensity pulse loads on actual structures. Figure 1; 
refere ices 10 (Russian). 


Threshold Diameter and Completeness of 
Detonation Decomposition of Coarsely Dispersed 
Commercial Explosives 


957A1192A Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 3-13 


[Article by K.K. Shvedov and V.V. Patsyuk, Chemical 
Physics Institute, Chernogolovka, Russian Academy 
of Sciences; manuscript received 13 Oct 94, UDC 
662.215.1]} 


([FBIS Abstract] A number of widely used commercial! 
coarse-grained individual explosives and mixed ammo- 
nium nitrate-based explosives were studied to determine 
the effect of the explosives’ particle size on their to- 
tal chemical decomposition and on the interaction of 
their components. Pairs of the following explosives were 
used: powdered trinitrotoluene [TNT]-granuloto!, 6ZhV 
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ammonite-79/21 grammonite, and a mixture of ammo- 
nium nitrate eid diesel fuel in a 94:6 ratio (igdanite). 
The explosive used within one and the same series of 
experiments had identical chemical or component com- 
positions; only their particle diameter (4) was varied. 
The experiments were performed with charges in card- 
board shells and steel tubes with diameters up to 200 
mm and wall thicknesses of 3-20 mm and with deep- 
hole charges with charge diameters [d] of 250-400 mm. 
In each experiment, two or three measurements of the 
detonation rate [D] at different distances from the ini- 
tiating explosive in the second half of the charge were 
generally taken. At least two or three experiments were 
performed for each diameter studied. Blocks of pressed 
TNT with a density [p] of 1 £9 g/cm? and 6 = 80 mm 
were used to initiate the explosions. The maximum error 
in determining D given a base of 100 mm did not exceed 
2 percent, and the variance from the mean amounted 
to 1 percent. Conclusions regarding the existence of a 
threshold detonation diameter were based on measure- 
ments of the dependence D(d) taken over a wide range 
of variations in charge diameter and on estimates. The 
existence of a threshold detonation diameter was shown 
for all of the explosives studied. The experiments estab- 
lished that in the case of the individual explosives, the 
ideal detonation rate does not depend on particle size. 
In the specific case of TNT with p = 0.9 g/cm’ in a 
charge with d = 145 mm in a steel tube, D = 4.62 km/ 
s. The same detonation rate was achieved with charges 
of granulotol with 4 = 1-5 mm and d = 80 mm in a 
steel shell and remained virtually constant as the charge 
diameter was increased to 145 mm and the thickness of 
the steel tube’s wall was increased from 3 to 20 mm. 
In the case of mixed explosives based on ammonium 
nitrate, on the other hand, the maximum detonation rate 
was indeed determined by the particle size: The more 
disperse the explosive (which is to say, the smaller its 
particle size), the greater ‘he detonation rate. The re- 
lationship between the detonation rate and dispersion 
(particle diameter) was concluded to be responsible for 
the incompleteness of the decomposition of mixed ex- 
plosives as a result of problems with diffusion given 
the high pressures behind the shock front. The discrep- 
ancy between the results yielded by actual experiments 
and those yielded by the existing theory of detonation 
and methods for computing detonation parameters and 
the effects of explosions was thus concluded to be due 
to the fact that the characteristic equations for detona- 
tion products and the physical principles on which the 
computation methods are based do not accurately reflect 
the physicochemical processes occurring in the reaction 
zone of the detonation wave. Figures 4, tables 3; refer- 
ences 24: 12 Russian, 12 Western. 
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Effect of the Shape of the Initial Shock Wave Pulse 
on the Decomposition of Polycrystalline Explosives 


957A1192B Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 14-26 


[Article by $.G. Andreyev, Moscow State Technical 
University imeni N.E. Bauman; manuscript received 29 
Sep 94; UDC 534.222.2+662.215.1] 


[FBIS Abstract] Published data were used to demon- 
strate the possibility of using principles and data regard- 
ing the viscoplastic mechanism of the formation and 
laminar decomposition of quasi-spherical seats of de- 
composition to predict the behavior of structurally het- 
erogeneous explosives given differently shaped initial 
shock wave pulses. Next, a series of original experi- 
ments were conducted to show the effect of the com- 
position and structure of explosives and the shape of 
the dynamic initial pulse on the unique gas-dynamic 
and electrical manifestations of decomposition of the 
elementary cell of the explosive charge in weak shock 
waves. In the experiments, an explosive that had been 
compacted into a table with a diameter of 40 mm and 
thickness of approximately 0.4 mm was sandwiched be- 
tween two electric plates (the bottom one had a diam- 
eter of 50 mm and the top one had a diameter of 18 
mm). The explosives and cover plates were pressed be- 
tween two blocks. The top block (which was 15 mm 
thick and made of fluoroplastic) had a manganese pri- 
mary detector attached to its bottom edge. The lower 
block served to generate the initial pulse after being 
struck by a 50-mm-diameter, 30-mm-thick aluminum 
striker that was propelled through the air gap by the 
detonation products of a charge of a low-density mix- 
ture of hexogen and mipor. The device generating the 
initial pulse was mounted 25 mm from the end of the 
100-mm-long accelerating tube. Single- and multiple- 
step pulses were generated. Tablets of pentaerythritol 
tetranitrate [PETN] and BTF were used. The experi- 
ments confirmed that the shape of the initial pulse does 
indeed have an effect on decomposition of the explosive 
PETN. The experiments also demonstrated the presence 
of effects, namely, inhibition of decomposition as a re- 
sult of an increase in external pressure on the reacting 
cell and acceleration of decomposition as a result of lo- 
cal liquefaction, that are difficult to explain in terms of 
conventional models of focal reactions. Specifically, it 
was discovered that in those charges studied with a den- 
sity equal to approximately 90 percent of the maximum 
possible density as a result of compression by a shock 
wave with a pressure of 1-3 GPa, two groups of foci 
of decomposition form. The first group of foci is sur- 
rounded by condensed matter in which discontinuities 
capable of significantly changing their gas permeability 
have not developed during the course of their expan- 
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sion. During the process of their development, the foci 
of the second group are surrounded by a gas-permeable 
or potentially gas-permeable medium and are located 
primarily on the intergrain boundaries and courses of 
cracks in the grains of explosive that were formed dur- 
ing the process of making the charge or else during the 
course of generation and development of the reaction 
centers. The said types of foci were said to be similar 
to foci that have been observed in direct optical exami- 
nations of polished sections of TNT charges after shock 
wave loading at pressures of 0.4-1.0 GPa. The mecha- 
nism of the interaction of the reaction products with the 
explosive in the vicinity of discontinuities in the ma- 
trix material in the initial stage of development of focal 
decomposition was discussed in detail. Figures 7; refer- 
ences 21: 16 Russian, 5 Western. 


Shock Wave Sensitivity of Explosives Irradiated 
With Soft X-Radiation 


957A1192C Moscow KHIMICHESKAYA FI 'KA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp. 28 


[Article by V.M. Belskiy, V.N. German, and IP. 
Khabarov, All-Russian Scientific Research Institute of 
Experimental Physics, Arzamas-16; manuscript received 
13 Oct 94; UDC 662.215.4] 


[FBIS Abstract] The "gap-test" method was used to de- 
termine the effect of irradiation with soft x-rays on the 
shock wave sensitivity of the hexogen-based explosive 
TG 30/70. Pieces of the explosive with a diameter of 
40 mm and height of 20 mm were bombarded with soft 
x-rays from URS-4 and ARS-4 units at dose rates of 3 
and 100 roentgens [R)/min, respectively. The dose of 
radiation directly absorbed by the material was not de- 
termined; rather, the degree of irradiation was judged 
on the basis of exposure. The exposure dose was varied 
over time, with the density of the study specimens re- 
maining the same before and after irradiation. All spec- 
imens were detonated within 24 hours of having been 
irradiated. As the exposure dose was increased, the deto- 
nation delay time initially decreased and then increased. 
Such shock wave sensitivity was said to be analogous 
to previously reported behavior and was explained in 
terms of the hypothesis of the development of foci of 
chemical reaction as a result of annihilation of mobile 
dislocations durine shock compression of the explosive. 
When the exposure dose used in the experiments was 
increased, both ‘he concentration of mobile dislocations 
in the grains of explosive and the concentration of dislo- 
cations of the foci of chemical reaction that develop as a 
result of annihilation increased. As a result, the detona- 
tion delay time decreased. After a specified threshold, 
however, the concentration of mobile dislocations de- 
creased as a result of blocking thereof. This in turn led 
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to a decrease in the concentration of foci, an increase 
in the detonation delay time, and the appearance of a 
trough on the curve plotted to represent the dependence 
of the detonation delay time of TG30/70 on the ex- 
posure dose. When the starting explosive already con- 
tained the maximum possible concentration of mobile 
dislocations, irradiation resulted solely in a decrease in 
the said concentration and, consequently, in an increase 
in the detonation delay time as compared with the start- 
ing detonation delay times. The experiments established 
that when all other conditions (charge density, amount 
of impurities, grain size, etc.) are equal, the degree of 
distortion of the crystalline lattice of the explosive itself 
is what determines an explosive’s shock wave sensitiv- 
ity. Figure 1; references 4: 3 Russian, | Western. 


Connection Between Explosives’ Critical and 
Maximum Diameters of Detonation and the 
Nonstationary Segment of Detonation With Shock 
Wave Sensitivity 


957A1192D Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 29-31 


[Article by I.P. Khabarov, V.M. Belskiy, and L.V. 
Fomicheva, All-Russian Scientific Research Institute of 
Experimental Physics, Arzamas-16, manuscript received 
13 Oct 94; UDC 662.215.4] 


([FBIS Abstract] The connection between the critical and 
maximum detonation diameters of an explosive and the 
nonstationary segment of detonation with shock wave 
sensitivity was examined in a series of experiments 
performed with specimens of compacted TG 30/70 
with an initial density of 1.720 +/- 0.003 g/cm’. The 
explosive’s critical thickness was determined based 
on attenuation of the detonation wave in truncated 
pyramids with an angle measuring of 8 to 2° at 
the vertex and a length of 120 mm. The explosive 
was detonated by an auxiliary detonator measuring 
10 mm in both diameter and height and made of 
phlegmatized hexogen. Cessation of detonation was 
judged on the basis of the trace left on the substrate’s 
surface as a result of the detonation wave's brisant 
effect. The error in measuring the critical diameter of 
detonation did not exceed 0.1 mm. The explosive's 
shock wave sensitivity was evaluated on the basis 
of detonation delay time as determined by the “gap- 
test” method. The maximum detonation diameter was 
determined from the dependences of the detonation 
rate on the charge diameter and on the length of that 
segment of development of the detonation. Detonation 
rate and lenght of the segment of development of 
the detonation were in turn detemrined by ionization 
detectors (in cases of a measurement base of 100-170, 
the measurement error did not exceed 0.20 percen 
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and in cases of a measurement base of 10 mm, the 
error in determining the detonation rate did not exceed 
1 percent). The experiments established that as the 
shock wave sensitivity of the explosive TG 30/70 
is decreased, its critical and maximum diameters of 
detonation and the length of the process of development 
of the detonation increase. In charges with a fixed 
diameter, the detonation rate was found to be a linear 
function of the detonation delay time. The slope of 
the line plotted to represent the said linear dependence 
became greater as the charge diameter became smaller. 
The observed decrease in detonation rate was attributed 
to an increase in the critical and maximum diameters 
of detonation. It was recommended that the aforesaid 
results be taken into consideration when the processes 
of excitation and development of a detonation process 
are studied and when initiation systems are developed. 
Figures 4; references 2: 1 Russian, | Western. 


Theoretical Substantiation of a Phenomenologic 
Model of the Shock Wave Sensitivity of a 
Heterogeneous Explosive Based on TATB With 
Consideration for One- and Two-Time Shock Wave 
Loading, Including With Intermediate Relief 


957A1192E Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 32-39 


[Article by V.G. Morozov, I.1. Karpenke, S.Ye. 
Kuratov, S.S. Sokolov, B.N. Shamrayev, and L.V. 
Dmitriyeva, All-Russian Scientific Research Institute 
of Experimental Physics, Arzamas-16; manuscript 
received | Nov 94; UDC 532.593] 


{[FBIS Abstract] A phenomologic model of the shock 
wave sensitivity of a heterogeneous explosive based 
on TATB is proposed that gives consideration to one- 
and two-time shock wave loading, including modes of 
shock wave loading with an intermediate stage of un- 
loading. Tie model is based on the existing concepts 
that “hot spots” are generated and develop on pores in 
heterogeneous explosives when a shock wave passes 
through them and that the said hot spots then transform 
the shock wave into a detonation wave. The innova- 
tion of the newly proposed model lies in the fact that 
it includes macroscopic equations of the kinetics of de- 
composition of an explosive. Specifically, the evolution 
of a hot spot is analyzed based on estimates and ap- 
proximate solutions describing the course of reactions 
in a “spot” with consideration for heat conduction and 
gas dynamics. The overall form of the rate of chem- 
ical reactions is varied without detailing their micro- 
scopic mechanism. The approach on which the newly 
proposed model is based is said to make it possible to 
obtain a type of macroscopic equation and, at the same 
time, draw conclusions regarding the characteristic fea- 
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tures of microscopic processes occurring within a hot 
spot. Specifically, the new model makes it possible to 
conclude the following: |) a strong increase in combus- 
tion and heat conduction-based energy transfer occurs 
in the form of a single wave and is linked by the com- 
mon mechanism of a change in molecular and electron 
structure of the material and 2) heat conduction plays a 
decisive role in the detonation process and depends on 
the pressure in the stage during which the shock wave 
becomes a detonation wave. It was hypothesized that 
a similar mechanism occurs during layer-by-layer com- 
bustion and dictates the pressure dependence of veloc- 
ity of the combustion front. The newly derived kinetics 
equations have been inserted in programs for numeri- 
cal calculation of the gas dynamics of the process of 
explosion of an explosive. Work to further develop the 
model of the evolution of hot spots is continuing in the 
following directions: investigation of nonequilibrium re- 
actions, especially at the initiation point and near the hot 
spot’s surface where the deformation rate is at its peak, 
and study of the evolution of hut spots on the surface 
of microgranules. Figures 4; reference 1 (Russian). 


Shock Wave Synthesis of Stannous Sulfide 


957A1192F Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 40-48 


[Article by §.S. Nabatov, G.Ye. Ivanchikhina, A.V. 
Kolesnikov, A.V. Kulbachevskiy, and A.V. Lebedev, 
Chemical Physics Institute, Chernogolovka, Russian 
Academy of Sciences: manuscript received 17 Oct 94; 
UDC 662.216.5] 


[FBIS Abstract) Published data regarding shock wave 
synthesis of stannous monosulfide [SnS)} were reviewed, 
and an original experiment was performed in which 
a special flat storage vial was used to take electrical 
measurements during the course of the reaction Sn + § 
— SnS in accordance with the shock wave synthesis 
technique. The publications reviewed were said to 
demonstrate that SnS may be synthesized in a shock- 
compressed reaction mixture. The original experiments 
were conducted in flat vials made of No. 3 steel that 
had a wall thickness of 5-10 mm on the side from 
which the shock wave entered. In the experiments, a 
Study specimen measuring 20 mm in diameter and 2-3 
mm in height was placed between the inner surface of 
the vial’s front wall and a massive internal electrode 
(diameter, 20 mm) that was isolated from the vial’s 
body by an insert and disk that had been turned from 
lithographic stone. The resultant structure was a plane- 
parallel electrical cell with the vial’s body serving as 
the grounded electrode. Shock wave loading of the 
Study specimen was accomplished by using a 90-mm- 
diameter aluminum striker accelerated by the explosion 
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products. A 7-mm-thick striker generated a rectangular 
wave with an amplitude of 13 GPa in the vial wall. A 
5-mm-thick striker generated a loading pressure of 16 
GPa. In some cases, a plane detonation wave (diameter, 
100 mm) generator that was in direct contact with the 
vial was used and resulted in a pressure of 8 GPa. 
In the range of pressures examined, the vial made 
it possible to make reliable electrical measurements 
both during dynamic effects and during the course of 
some arbitrarily long time interval after arrival of the 
unloading wave. Several oscillographs with different 
sweep durations (from approximately 10 ys to 1 s) were 
used to continuously record the electrical parameters 
of the material inside the vial. The study specimens 
were prepared from a mixture of tin and sulfur powders 
(grain size, 5-20 jam) that were used in a stoichiometric 
ratio, mixed, and compacted to a density of 4.3 g/ 
cm’. When the vials were loaded at pressures of 13 
and 6 GPa, the retained products consisted virtually 
entirely of tin monosulfide. No reaction was initiated 
at a pressure of 8 GPa. After the published and original 
results were analyzed together, it was concluded that 
depending on the intensity of the dynamic effect, the 
reaction Sn + § —+ SnS may occur not just in a shock- 
compressed reaction mixture but also in an unloaded 
reaction mixture. A method of recording the reaction 
Sn + S — SnS behind the shock wave front based on 
the electrical response generated by a semiconductor 
compound during the shock compression was proposed. 
Figures 5; references 17: 14 Russian, 3 Western. 


Unloading Isentropes and Characteristic Equation 
of Molybdenum at High Energy Densities 


957A1192G Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 49-55 


{Aciscie by M.V. Zhernokletov, A.B. Medvedev, and 
G.V. Simakov, All-Russian Scientific Research Insti- 
tute of Experimental Physics, Arzamas-16; manuscript 
received 13 Oct 94; UDC 532.593+536.2] 


(FBIS Abstract] A study examined the behavior of dense 
molybdenum plasma in shock compression waves and 
its subsequent isentropic unloading. The states in the 
unloading waves were recorded by using the method 
of adiabatic expansion of the shock-compressed metal 
into media with lower dynamic impedance values and 
known shock adiabatic lines. An enclosure with a 
different rigidity was used so as to make it possible 
to follow the course of the expansion isentrope at the 
coordinates pressure (P and mass velocity (U) from 
states on the shock adiabatic line from the region 
representing pressures of hundreds of gigapascals all 
the way to the gas region, where pressures amount to 
tens of megapascals, and to determine the change in 
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specific volume and specific energy of the expanding 
metal based on differential Riemann relauonships. Two 
series of experiments were performed. The molybdenum 
specimens used in the first series had a density of 
4.43 g/cm’, which corresponds to a porosity of 2.3. In 
the second series of experiments, scientists studied the 
unloading isentropes of contuinuos (dense) specimens 
(with a normal density of 102 g/cm’) from three 
states on the shock adiabatic line with a porosity of 
1. The parameters of the isentropes’ initial states were 
determitied by the reflection method. The measurements 
taken during the experiments were plotted along the 
coordinates P-U. The resultant curve confirmed that the 
properties of molybdenum change continuously as it 
expands from a condensed state without any noticeable 
sharp bends that could be interpreted as specific plasma 
phase transitions. The only anomaly occurred on the 
unloading isentropes of the porous specimens when 
they entered the near-critical region, which resulted 
in an additional increase in the rate of expansion and 
permitted qualitative determination of the boundary of 
high-temperature boiling. The experiment results were 
compared with the results of calculations based on a 
characteristic equation that takes the effects of melting, 
vaporization, and ionization into account. The two sets 
of results were found to be in good agreement with one 
another. It was determined that during shock loading, 
molybdenum begins at a pressure of 403 GPa and ends 
at a pressure of 493 GPa. Figures 4, tables 2; references 
14: 7 Russian, 7 Western. 


Characteristic Equation of Liquid and Gaseous 
Cesium 


957A1192H Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 56-64 


[Article by A.B. Medvedev, All-Russian Scienufic Re- 
search Institute of Experimental Physics, .arzamas-16; 
manuscript received 13 Oct 94; UDC 536.71] 


{[FBIS Abstract] A modified Van der Waals model was 
used to construct a characteristic equation for liquid 
and gaseous cesium that reflects the presence of va- 
porizatuion and thermal ionization. The model consid- 
ers only the thermal ionization that is characteristic for 
gaseous systems but not ionization “by pressure,” and 
the metal-dielectric transition is ignored. The contribu- 
tion of conduction electrons to compressibility is effec- 
tively taken into consideration by the empirical covol- 
ume (incompressible volume). According to calculations 
based on the modified model, cesium’s temperature in- 
creases greatly when it is subjected to shock compres- 
sion. Specifically, model temperature values of 27,000- 
45,000 K are obtained for pressures of 30-50 GPa. A 
comparison of the results yielded by calculations based 
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on the proposed model and results obtained by actual 
experiments established that the agreement between the 
two sets of results is fairly good within a broad range 
of parameters extending from threefold degrees of com- 
pression to thousandfold degrees of expansion and tem- 
peratures of several tens of thousands of degrees. The 
biggest discrepancies between calculated and experi- 
mentally obtained liquid-vapor equilibrium curves oc- 
curred in the region adjacent to the critical point (where 
T is greater than or approximately equal to 0.9 times 
the critical temperature), at which point fluctuations not 
taken into consideration in Van der Waals-type models 
make a significant contribution to cesium’s thermody- 
namic properties. Figures 6; references 16: 13 Russian, 
3 Western. 


Adiabatic Line of Uranium Up to Pressures of 4 
TPa 


957411921 Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 65-67 


[Article by L.V. Altshuler, A.A. Bakanova, M.I. Brazh- 
nik, V.I. Zhuchikhin, $.B. Kormer, K.K. Krupnikov, and 
R.F. Trunin, All-Russian Scientific Research Institute of 
Experimenta! Physics, Arzamas-16; manuscript received 
13 Oct 94; UDC 532.593] 


[FBIS Abstract) Research on uranium’s shock com- 
pressibility began in 1947 at Russia’s Arzamas Nuclear 
Center and continued through the 1950s and 1960s. In 
1966, uranium’s shock compressibility under pressures 
of 4.0 terapascals was studied in an underground explo- 
sion at Novaya Zemlya. Uranium disks with a thickness 
of 40 mm were placed behind a lead screen and sub- 
jected to an explosion with a shock wave having a peak 
velocity of D,, mF 19.1 km/s. Wave velocity measure- 
ments were taken as the shock wave passed through 
Pb-U samples. The characteristics of uranium’s shock 
compressibility were determined by plotting a pressure- 
mass velocity [P-U] diagram based on the known pair 
of values D, D,, and the shock adiabatic line of lead 
(which, in the range of values studied, corresponds to a 
value of D = 3.19 + 1.167U, with U representing mass 
velocity and assuming that the initial density of the lead 
equaled p, = 11.34 g/cm’). The results were compared 
with results obtained in three Western laboratories and 
published in 1962 and 1980. The results obtained at No- 
vaya Zemlya were said to be in fairly good agreement 
with those obtained by Western scientists and reported 
in a 1980 publication. The consistency of the results 
obtained at Novaya Zemlya, which are of an absolute 
nature, and the Western results, which are of a com- 
parative nature, was said to be confirmation of the cor- 
rectness of both approaches. The Russian results were 
said to reflect rather large measurement errors in the 
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range of high parameters (D imF 10.5S-13.0 km/s) that 
were attributed to limitations of the measurement equip- 
ment used when the measurements were actually taken 
(1951-1953), including excessive heating of the strikers 
and unsatisfactory (by today’s standards) symmetry of 
the flight of the striker to the target. It was emphasized, 
however, that the average D-U line is sufficiently reli- 
able and has been validated by control measurements 
taken with modern equipment (at D ImF 10.5 km/s). 
Figure 1; references 3 (Western). 


Shock Compressibility of Iron, Aluminum, and 
Titanium at Terapascal Pressures 


957A1192J Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 97-99 


[Article by R.R. Trunin, N.V. Panov, and AB. 
Medvedev, All-Russian Scientific Research Institute 
of Experimental Physics, Arzamas-16; manuscript 
received 13 Oct 94; UDC 532.593} 


([FBIS Abstract} For the first ume, data have been ob- 
tained regarding the shock compressibility of aluminum 
and tantalum under pressures of 0.689 and 1.66 terapas- 
cals [TPa] on a laboratory measuring unit. The measure- 
ments were taken with iron, which has a shock compres- 
sion pressure of 1.26 TPa, serving as a standard. The 
unit used to conduct the experiments is a cumulative de- 
vice in which the products of a strong explosion accel- 
erate a specially shaped steel striker-liner to a velocity 
of 18.5 km/s with the symmetry of the shock wave cre- 
ated in the study specimens upon collision with the liner 
being no worse than 10° s and with the initial heating 
of the striker by the shock wave being approximately 
equal to the effect that the planar-wave detonation front 
has on it (shock wave amplitude, 45 GPa). Several tests 
measuring the shock compressibility of samples of Fe, 
Al, and Ta with initial densities of 7.85, 2.71, and 16.65 
g/cm’, respectively, were conducted. In the tests, sam- 
ples with a diameter of 10 mm and thickness of 4 mm 
were placed behind a 3.5-mm iron screen. The follow- 
ing average wave velocities were established for Fe, Al, 
and Ta, respectively (km/s): 17.35 +/- 0.3, 20.9 +/- 0.4, 
and 13.5 +/- 0.9. All of the said values include small ( 
2 percent) positive corrections to allow for the incon- 
stancy of velocity in the samples. The following values 
were calculated for the mass velocity [U], pressure [P], 
density (], and degree of compression (0} of aluminum: 
U = 12.20 km/s; P = 0.689 TPa; p = 6.3 g/cm’: and o = 
2.41. The corresponding values calculated for tantalum 
were as follows: U = 7.48 km/s; P = 1.656 TPa; p = 
38.05 g/cm’; and 7 = 2.29. The aforesaid results were 
found io be in satisfactory agreement with data obtained 
at lower pressures, data obtained in experiments using 
strong shock waves in underground explosions (in the 
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case of Al), and shock adiabatic lines calculated on the 
basis of a model characteristic equation. Figure 1; ref- 
erences 11: 10 Russian, 1 Western. 


Effect of a High Pulse Load on Massive Metal Disks 


957A1192K Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 112-116 


[Article by V.I. Zeldovich, N.P. Purygin, B.V. Litvinov, 
1.V. Khomskaya, G.A. Sobyanina, N.D. Matushkin, and 
V.V. Medvenko, All-Russian Scientific Research Insti- 
tute of Technical Physics, Snezhinsk (Chelyabinsk-70); 
manuscript received 4 Nov 94; UDC 539.89+620.18] 


(FBIS Abstract] Massive metal disks were subjected to 
shock waves that exerted a pressure of approximately 
100 GPa on their centers. The deformation, crack 
formation, and microstructure of the massive disks were 
examined after they had been subjected to convergent 
shock waves in an axisymmetric pulse loading mode. 
Three materials, i.e., type St3 steel, VTZ-1 titanium, 
and D16 duralumin, were studied. Each material was 
fashioned into a continuous cylindrical disk with a 
diameter of 60 mm and height of 30 mm. The pressures 
along the axis of the steel, titanium, and duralumin disks 
when the detonation wave reached their face amounted 
to 560, 480, and 420 GPa, respectively, and the pressure 
in their centers after collision with the shock wave 
equaled 130, 110, and 100 GPa, respectively. When 
subjected to shock wave loading with a pressure of 
approximately 100 GPa, the duralumin disks proved 
much more ductile than the steel and titanium disks. The 
deformation that developed in the titanium disks upon 
collision with the shock wave was much more localized 
in the center of the disks than was the case with the 
steel disk. The steel disks contained fewer “surface” 
cracks than the duralumin or titanium disks; however, 
the steel disks contained a round crack located in their 
central piane that was not present in the other two types 
of disks. In the steel disks, cracks formed as a result 
of tensile stresses that developed upon collision with 
the shock wave. In the titanium disks, cracks formed 
as a result of shear stresses and were a continuation of 
adiabatic shear bands. A link between crack location and 
the direction of propagation of the direct and reflected 
shock waves was discovered in the case of the steel 
disks. Figures 4, table 1; references 8: 7 Russian, | 
Western. 


Phase Transition in Quartzite and Scale Effect 


957A1192L Moscow KRIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 122-125 


[Article by V.G. Vildanov, M.M. Gorshkov, V.M. Slo- 
bodenyukov, and P.N. Senichev, All-Russian Scien- 
tific Research Institute of Technical Physics, Snezhinsk 
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(Chelyabinsk-70); manuscript received 1 Nov 94; UDC 
536.42+539.89] 


[FBIS Abstract] The structure of shock wave fronts de- 
veloping in virtually pure quartzite (96 percent SiO, 
with a starting density of 2.65 g/cm’) was examined in 
a series of experiments during which wave velocity and 
free surface velocity were measured under conditions of 
a strong underground explosion. During the explosion, 
several tens of kilotons of energy were released so as 
to make it possible to take measurements in large mea- 
surement bases with P 2 21.6 GPa. Electromechanical 
sensors were used to record the shock wave’s motion 
throughout the quartzite at the normalized distances | 
<r < 3.34 m/kg'” The measurement base (distance be- 
tween adjacent levels) amounted to Ar = 0.1 m/kt'” and 
was determined with a precision of at least 0.5 percent. 
Two sensors were mounted at each distance r, from the 
center of the explosion (level) for the sake of redun- 
dancy. The error in measuring time intervals did not 
exceed 0.1 percent, which means that the error in de- 
termining wave velocity equaled approximately 0.5 per- 
cent. Polished walls with a nonplanarity not exceeding 
+/-0.0002 m/kg'” were erected perpendicular to the di- 
rection of the shock wave’s propagation to determine 
the hodograph of the free surface H(t) at distances of r, 
= 3.01 m/kt’® and r, = 3.34 m/kt’” (the two distances 
corresponding to the region of the mixture of quartz 
phases) directly in the quartzite. The error in determin- 
ing the hodograph A(t) amounted to 0.5 percent. The 
measurements taken at the first free surface were used 
to establish the complex profile of the shock wave front 
consisting of a jump with an amplitude of approximately 
22 GPs, compression wave, and subsequent jump. At- 
tenuation of the shock wave in quartzite was determined 
to be described by the function D(r) = 20.47 (1/r)'?°, 
with D expressed in km/s and r expressed in m/kt'” 
(with a mean-square approximation error of § = 0.18 
km/s). When r 2 2.77 m/kt'’® D became constant and 
equaled 5.72 km/s. The existence of a second shock 
wave was detected when W(t) was measured at a dis- 
tance of 0.24 m/kt'” from the point of the sharp bend in 
the dependence D(r) and with a recording time of 14 us/ 
kt'®. The resultant wave configuration was determined 
to be nonstationary: The velocity of the second shock 
wave was significantly less than the velocity of the first 
(which remained constant at D = 5.72 km/s). Figures 2; 
references 8: 4 Russian, 4 Western. 
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Methodological Features of the Study of Dynamic 
Compressibility of Soft Ground in the Range of 
Pressures From 0.05 to 1.5 GPa 


957A1192M Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 126-132 


[Article by A.M. Bragov, V.P. Gandurin, G.M. Gn- 
shevskiy, and A.K. Lomunov, Mechanics Scientific Re- 
search Institute, Nizhniy Novgorod University; manu- 
script received 17 Oct 94; UDC 534.222.2:624.131] 


[FBIS Abstract] Two similar methods of testing soft 
ground under laboratory conditions are described that 
make it possible to study dynamic compressibility in a 
previously unresearchable range of pressures, namely, 
pressures in the range from 0.05 to 1.5 GPa. At 
pressures up to 300 MPa, a modification of the Kolskiy 
method of using a Hopkinson split bar with passive 
limitation of radial deformation is used. Specifically, the 
proposed method makes it possible to construct dynamic 
diagrams of soft soils under conditions of uniaxial 
deformation and also to determine their cubic strain 
curve and the pressure dependence of their shearing 
strengt. The main difference between the setup of the 
proposed method and the conventional Hopkinson split 
bar method is that in the proposed method the soil 
sample is placed in a rigid steel casing that prevents 
it from undergoing any radial deformation. At pressures 
above 300 MPa and deformation rates of more than 10*/ 
S, a plane wave impact experiment is used that entails 
a gas gun with a caliber of 57 mm. The process of 
using the proposed methods is illustrated by way of 
the example of the construction of dynamic diagrams 
of the compressibility of real soils (clay and sand) 
and plasticine. Four or five tests of each sample were 
conducted under identical conditions, and the results 
were averaged. The deformation curves plotted for 
plasticine were clearly nonlinear in their loading leg. 
In the range of deformation rates studied, the initial 
segments of the deformation curves plotted for all four 
plasticine samples were nearly identical. The sections of 
the four samples’ curves in the legs reflecting unloading 
were all different from one another and were determined 
to be largely dependent on the deformation rate and 
levels of stress or strain achieved during the individual 
experiment. A comparison of the deformation diagrams 
plotted for plasticine and clay soil with a density of 2 
kg/cm’, moisture content of 24 percent, and porosity 
of 43 percent revealed that when subjected to active 
loading, plasticine’s behavior mimics that of actual 
clay soil rather closely. A comparison of the results 
of experiments in which the clay soil was subjected 
to dynamic and static loading established that clay's 
compressibility is significantly reduced under dynamic 
loading. The plane wave experiments established a 
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linear relationship between shock wave velocity and 
mass velocity of the soil and plasticine particles. Figures 
8, table 1; references 9: 5 Russian, 4 Western. 


Refined Model of Plastic Deformation of a Porous 
Medium 


957A1i92N Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 136-142 


{Article by A.I]. Sadyrin, Mechanics Scienufic Research 
Institute, Nizhniy Novgorod University; manuscript re- 
ceived 19 Oct 94; UDC 539.3:621.7] 


[FBIS Abstract] The so-called Herrman-Carrol-Holt P-a 
model is generally used to analyze the results of experi- 
ments studying the shock wave loading of many porous 
media and to calculate the deformation of porous ma- 
terials in structural elements. The model has been re- 
fined for the case of moderately porous media (a > 
G.5) with loads not significantly exceeding the thresh- 
old of complete closing of the pores. The dependence 
of pressure in the porous medium on the specific vol- 
ume, ie., P(v), is represented in parametric form by 
introducing the parameter a that characterizes the mate- 
rial’s degree of porosity through two functions: r. = flv, 
a) and a = g(P,). The newly derived relationships are 
compared with existing expressions for calculating plas- 
tic deformation during isotropic compression. The new 
relationships are generalized within the framework of 
an associated theory of plasticity to give consideration 
io the deviator components of stresses. A comparison 
of the calculated density dependence of the elasticity 
modulus of AL,O,-based porous ceramic that was per- 
formed in accordance with the proposed model and ex- 
periment results reported elsewhere confirmed that the 
proposed model yie.ds results consistent with the re- 
sults of experiments. In addition, a comparison of the 
calculated changes in the intensity of stresses during 
uniaxial deformation and a uniaxial stressed state dur- 
ing the compression of porous aluminum (a,, = 0.78) 
and the results of actual experiments established that the 
proposed refined relationships describe selected quali- 
tative and quantitative characteristics of the experimen- 
tally observed plastic deformation of porous media with 
satisfactory precision. Figures 3; references 20: 9 Rus- 
sian, 11 Western. 


Excitation of a Turbulent Flow by 
Richtmayer-Meshkov and Rayleigh-Taylor 
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Instability During a Cylindrical Convergent 
Explosion 

957A11920 Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 150-158 


[Article by M. Legran and M. Nikolayzo, Atomic En- 
ergy Commissariat, Research Center, France; manu- 
script received 29 Aug 94; UDC 532.593) 


(FBIS Abstract] The excitation of a turbulent flow by 
Richtmayer-Meshkov and Rayleigh-Taylor instability 
during a cylindric convergent explosion was examined 
in an experiment. A convergent explosion wave was 
initiated by a cylindrical detonation wave propagated 
along a cylindrical explosive charge (outer diameter, 
200 mm). At the outset of the experiment, a defect 
was created on the inner surface of the tin shell 
(inner diameter, 100 mm; wall thickness, 4 mm) that 
represented a combination of sinusoidal waves. The 
plan of the experiment was to have the Rayleigh-Taylor 
instability develop over a lengthy period as a result of 
the rarefaction wave moving behind the detonation wave 
and counterpressure developing in the low-density inner 
medium. A one-dimensional hydrodynamic code was 
used to calculate the time dependence of the velocity 
of the interface between the tin and silicon in the 
central part of the cylinder. The experiment results 
were similar to the results of numerical calculations 
based on a two-dimensional Euler code; however, the 
flow was hardly turbulent inasmuch as the area near 
the interface was not a mixing zone. Attempts to 
modify the spectrum of initial defects in the tin were 
moderately successful: The experiment results turned 
out to be very close to the calculated results; however, 
the nature of the flow (which could have just as likely 
been laminar as turbulent) could not be determined 
with certainty. Thus the experiment demonstrated the 
development of different disturbances only during the 
course of a preiurbulent phase. Next, the criteria of the 
development of a turbulent layer near an interface were 
considered. It was demonstrated that turbulence arises 
dunng the initial disturbances of very high modes of 
instability. A generalized Youngs-Read model was used 
in a plane geometric approximation so as to make it 
possible to compare calculations in a plane geometry 
with the results of an experiment reflecting a cylindrical 
system. The increases in mixing zone width over 
ume calculated by numerical modeling (with the two- 
dimensional Euler code) and by using the generalized 
Youngs-Read model turned out to be very close to one 
another. When the Youngs-Read model was applied to 
the case of a cylindrical configuration, the mixing zone 
width increased much faster than in the case of a plane 
configuration. Figures 8; references 7 (Western). 
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Electrical Response of Flexible Piezoelectric 
Composite Films to High-Speed Collision With a 
Spherical Body 

957A1192P Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 159-165 


[Article by V.M. Shunin, $.S. Nabatov, L.A. Smetan- 
ina, and V.V. Yakushev, Chemical Physics Institute, 
Chernogolovka, Russian Academy of Sciences; manu- 
script received 1 Nov 94; UDC 539.89+637.226+678) 


([FBIS Abstract] The electrical response of piezoelectric 
composite films to high-speed collision with a spherical 
body was examined in a series of experiments invoiv- 
ing 4- to 10-nim-diameter steel spheres accelerated to 
velocities of 0.7-2.0 km/s. The pie7electric composite 
films were heat-resistant films of organosilicon binder 
(SKTN organosilicon rubber with modifying additives) 
filled with finely disperse TsTS-19 ceramic powder. The 
ratio of components was suvh that the ceramic powder 
constituted 35 percent of the mixture’s volume. After 
the mixture of binder and powder had been premixed, it 
was processed in an ultrasound disperser and applied to 
7-ym-thick aluminum foil. The polymerization process 
was completed within 48 hours. The solvent residues 
were removed by vacuum drying. The films studied had 
an area of 150 x 300 mm’ and thickness of approx- 
imately 200-400 jm. The finished films had an aver- 
age density of 3.70 g/cm’. A calculation of the density 
of the two-component mixture based on SKTN binder 
with modifying additives filled with TsTS-19 piezoce- 
ramic in a concentration of 35 percent yielded a value of 
3.95 g/cm?, which meant that the films studied could be 
considered three-component films (with air inclusions 
constituting the third phase). The films were polarized 
in an electric field with an intensity of 5-8 kV/mm and 
then heated to 120-150°C, held at that temperature for 
at least 30 minutes, and then cooled to room temper- 
ature. The films were then subjected to artificial aging 
at 150°C for 1 hour. A setup containing the study film 
(which functions as a piezoelectric composite sensor), 
two electrodes, coaxial cable, optical velocimeter, spher- 
ical striker, and load impedance was used to record the 
films’ response to the high-speed collision. The elec- 
trical voltages generated by the piezoelectric composite 
films were on the order of tens of volts at a resistance 
of 50 ohms and proved to have the best time resolution 
of all known collision sensors (approximately 1 us). 
Within the range of velocities studied, the amplitude of 
uhe electrical response was linked to the nature of the 
effect (i.e., to the velocity and area of the sphere’s cross 
section). The experiments thus confirmed the possibil- 
ity of using piezoelectric composite films as the active 
medium of wide-format reusable collision sensors. Fig- 
ures 4, table 1; references 9: 6 Russian, 3 Western 
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Stability of a Shock Wave Generated by a Point 
Explosion in Gas With an Exponential Dependence 
of Density on Radius 

957A1192Q Moscow KHIMICHESKAYA FIZIKA 

in Russian Vol 14 No 2-3, Feb-Mar 95 pp 166-170 


[Article by V.M. Ktitorov, All-Russian Scientific Re- 
search Institute of Experimental Physics, Arzamas-16; 
manuscript received 8 Aug 94; UDC 532.593) 


(FBIS Abstract] The stability of a shock wave generated 
by a point explosion in gas with a density that is an ex- 
ponential function of its radius was analyzed in a linear 
approximation. Specifically, stability was analyzed with 
respect to linear nonradial perturbations located along 
spherical harmonics Y_. The regularly varying func- 
tions considered described asymptotic conditions of the 
development of perturbation corresponding to time pe- 
riods significantly longer than the time periods required 
for a sound wave to pass through the entire region of 
the explosion. This meant that the perturbations are in- 
dependent of the initial conditions and also that the per- 
turbations in the parameters at the wave front are expo- 
nentially dependent on time. The analysis was begun by 
deriving linear equations for the perturbations by using 
a linearization of hydrodynamic equations written in La- 
grangian form. The problem formulated was shown to 
have an analytic solution in the case where the adiabatic 
exponent () is close to 1. It was further shown that in 
the case where 7 = 1, the eigenvalues A(n) were shown 
to satisfy a simple equation, whereas in the general case 
(y > 1), the problem must be solved numerically. The 
critical value of the adiabatic exponent 7, was shown to 
be a weak function of &. For values of k equal to -2, -1, 
0, 1, and 1, 7, was found to equal 1.12, 1.17, 1.2, 1.21, 
and 1.22, respectively. The first harmonic was examined 
separately. It was shown to have two eigenvalues that 
are both real and negative. It was stated that the larger 
of the two is independent of 7, can be found by using 
the formula A,(1) = -1.5 - &/2, and corresponds to the 
linear propagation of the explosion wave with conser- 
vation of the quantity of motion. Figures 3; references 
5: 2 Russian, 3 Western. 


Perturbations of a Stationary Shock Wave 
Propagated Along a Channel 


957A1192R Moscow KHIMICHESKAYA FIZIKA 
in Russian Vol 14 No 2-3, Feb-Mar 95 pp 171-173 


[Article by V.M. Ktitorov, All-Russian Scientific Re- 
search Institute of tal Physics, Arzamas-16; 
manuscript received 8 Aug 94; UDC 532.593) 


(FBIS Abstract) The propagation of a shock wave in 
channels whose inner walls have dynamic characteris- 
tics (initial density, characteristic equation, etc.) differ- 
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ent from those of their outer materials was analyzed. The 
analysis focused on the general case with arbitrary char- 
acteristic equations of both materials. It was assumed 
that a plane stationary wave was propagated along a 
channel (parallel to the channel walls) as a result of 
pressure applied by a hot piston. Two special cases were 
considered in detail, i.¢., the case where both materials 
were ideal gases with densities very similar to one an- 
other and the case of a cylindrical nonuniformly filled 
channel in a homogeneous external medium. Reference 
1 (Russian). 


Generation of Compression Waves in Dielectrics by 
Pulsed Electron Beam 


957A1202A Moscow KHIMICHESKAYA FIZIKA 
in Russian Jan 95 Vol 14 No 1, pp 100-107 


[Article by V.V. Milyavskiy and V.A. Skvortsov, United 
High Temperatures Research Institute, Russian Acad- 
emy of Sciences, Moscow; manuscript received 05 Sep 
94; first paragraph is KHIMICHESKAYA FIZIKA ab- 
stract; UDC 539.2;539.12.04) 


(FBIS Translated Text) A mathematical model has been 
constructed and used to perform a numerical study of the 
process of the interaction of an intensive relativistic elec- 
tron beam with a solid high-molecular-weight dielectric. 
The model, which is based on a system of equations of 
a continuous (elastoplastic) medium, equations of kinet- 
ics and charge carriers, and Poisson equations, makes 
it possible to describe shock wave disturbances and the 
evolution of an electrical field in an irradiated object 
with consideration for the accumulation and relaxation 
of a space charge. 


The question of the degree and nature of the effect of 
electrical fields that develop in a test object as a result of 
accumulation of a space charge affects the dynamics of 
shock wave disturbances and becomes important during 
interpretation of the results of experiments in which 
relativistic electron beams are used to study the dynamic 
loading of solid dielectrics. This article is devoted 
to studying the aforesaid aspect of the interaction of 
relativistic electron beams with a dielectric. 


Many interconnected physical processes are initiated 
when a relativistic electron beam interacts with matter: 
shock waves afe generated; mechanical breakdown of 
the matter, melting, ionization, and destruction of the 
molecules occur; free radicals form; a bulk charge 
accumulates in the object and electrical currents are 
induced; variable electrical and magnetic fields develop; 
and electrical breakdown is possible. 

The following three components of the problem of 
describing the interaction of relativistic electron beams 
with matter may be identified: 





[ie ear oor heed Coping wl Sm ——— 
| is prohibited without permission of the copyright owners. 











FBIS-UST-95-052 
19 December 1995 


— Calculating the distribution of the power density of 
the energy released by the relativistic electron beam 
throughout the depth of the absorber (Q.], the rate of 
volumetric injection of the electrons (the number of 
electrons in the beam that have been thermalized per 
unit volume per unit time (5S), and the rates of volumetric 
generation of electron-hole pairs [G); 


— modeling the processes of shock wave generation 
and mechanical breakdown of the irradiated object; 


— modeling processes of the accumulation and relax- 
ation of space charges, the evolution of electrical cur- 
rents and fields, and the development of electrical break- 
down. 


We will solve the problem of interaction of a relativistic 
electron beam with a dielectric in a one-dimensional 
approximation under the assumption that the object is 
covered on both sides by a thin grounded metallic foil. 
The thickness of the foil is assumed to be negligible in 
comparison to the thickness of the object; consequently, 
its effect on the process of absorption of electrons in the 
matter and hydrodynamic processes may be ignored. It 
is further assumed that the electrons fall onto the object 
in a direction normal to its surface. 


To calculate the functions Q.(h,t) and S(h,t) when 
an electrical field is present in the object, we will 
use a semiempirical method' based on the empirical 
method? and on the solution of the following equation 
describing the passage of electrons through matter that 
ignores scattering and that approximates continuous 
deceleration: 


Tae ~ 557) - EM, 


where T is the kinetic energy of the electron, h is the 
Lagrangian ccordinate, ¢ is the absolute value of the 
electron’s charge, p(h) is the density of the matter, p, is 
the initial density, E(A) is the intensity of the electrical 
field, and $,(7,) is the stopping power of the electrons 
in the matter.’ 

As designed, the method is characterized by the fol!ow- 
ing unique features: 

— in the case where E = 0, the results obtained are 
analogous to those reported elsewhere’, 

— it adequately describes the change in electron path- 
length in matter for all distributions E(h),; 

— it gives consideration to ‘he fact that part of the 
beam's energy becomes potential charge energy in the 
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electrical field (most of this energy is later released 
when the object is heated by conduction currents). 


Figure | presents the results obtained when this method 
was used to calculate the functions D(h) (Q.(A,1) = 
AgD(hVe, where J(t) is the density of the current of 








the incident electron beam).' 
Feil 
5 
P . 
0 
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Figure 1. Graphs of the functions of the specific 
energy release of electrons with an energy of 0.5 


MeV in polystyrene: 1, calculation based on 
Tabata’s method’; 2-4, calculations by the method 
described elsewhere! (2, E = 0.3; 3, E = 2 MV/cm 
with a constant decelerating field; 4, E = -0.5 MV/ 

cm with a constant accelerating field). (X-axis 
represents h in cm; Y-axis represents D in 
MeV cm’/g.) 


To model the processes of generation and propagation 
of compression and rarefaction waves in an object, we 
will use a system of equations of the mechanisms of 
a continuous compressible medium (a one-dimensional 
elastoplastic approximation )*: 


where V is the specific volume of the matter, U is 
the mass velocity of the matter, P is the hydrostatic 
pressure, s. is the xx-component of the stress deviator, 
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a, is the principal stress along the x-axis, ¢ is the specific 
internal energy, Q,(h) and Q,(h) are the specific power 
of the energy release as a result of deceleration of 
the beam's electrons in the matter and passage of the 
electrical current, and y is the shear modulus. 


Equation system (1) is supplemented with a wide-range 
equation of state for the matter, i.c., P = P(V, «),° and 
the condition of Mises fluidity.° 


The processes of accumulation and relaxation of a space 
charge in high-molecular-weight dielectrics is described 
on the basis of the band theory of a solid body. 
Most high-molecular-weight compounds have no long- 
range order in their arrangement of atoms; therefore, 
some of the concepts entailed in the classic theory of 
crystalline bodies lose their meaning. We will present 
selected principles of the band theory that will be used 
henceforth with respect to amorphous materials. 


First, the chemical bonds of adjacent molecules make it 
impossible to disrupt short-range order. Consequently, 
the band structure of an amorphous dielectric is not very 
different from the structure of a crystal with a simi- 
lar short-range order. Because of numerous lattice de- 
fects, however, strictly localized electron states (traps) 
develop in the forbidden band. The values of the en- 
ergies «. and «, (on the energy scale, they respectively 
separate the band of localized states from the conduc- 
tion band and the valence band) are termed mobility 
thresholds. 


Second, the concept of an electron-hole passive band’ is 
important. A fast electron gives up its energy primarily 
as a result of ionization by creating electron-hole pairs. 
Gradually, as it losses energy, it also loses its ability 
to ionize lower atomic levels. When the energy of the 
electron relative to the bottom of the conduction band 
falls below some threshold energy (¢_) that is of the 
same order of magnitude as the width of the forbidden 
band (¢,), it loses its ionization capability altogether. 
The same be made said for holes.'’ States of electrons 
(holes) with an energy in the interval «, belong to the 
electron-hole passive band; «, + ¢, (€, - €,; €,). As 
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a rule, the width of the valence band is less than te 
(¢, = €,), and it is entirely within the electron-hole 
passive band. After having entered the electron-hole 
passive band, an electron sharply reduces its pace of 
relaxation from an energy standpoint (because it loses 
energy in smaller portions) and is thermalized within 
a period on the order of 10°'* seconds.’ (Footnote: 
An electron may tentatively be considered thermalized 
when it has entered a narrow electrical interval with a 
width of approximately 0.01 eV (hv, = 0.1 eV being the 
characteristic energy of an optical phonon) around the 
bottom of the conduction band.) 


Third, it should be noted that electrons and holes are 
not generated in the electron-hole passive band given 
the relativistic electron beam parameters in which we 
are interested (see below).' In certain cases, this makes 
it possible to consider the nature of the motion and 
scatter of carriers from a phenomenologic standpoint 
as the motion of particles of an ideal gas.* 


Proceeding from the aforesaid model of a solid high- 
molecular-weight dielectric, we will in this article 
consider two possible mechanisms of conduction: 


— conduction resulting from electrons (holes) located 
in the interval « along the energy scale; ¢, , hu.(e, hv,; 
€,) (low-energy conduction); the dielectric’s intrinsic 
conductivity is low-energy'?; 


— conduction resulting from electrons (holes) located 
in the interval «. + hv, along the energy scale; «, + ¢,; 
(¢, “yp: €, “Av,) (high-energy conduction). High-energy 
charge carriers are only present in the study object while 
it is being irradiated.'’ 

The aforesaid concepts make it possible to write the 
following equation system to describe the processes of 
accumulation and relaxation of a space charge: 
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The following designations have been used: N and 
P are the total concentrations of electrons and holes 
in the study object; N, and P, N, and P,, and n 
and , designate the concentrations of free low-energy 
electrons and holes, electrons and holes fixed in traps, 
and high-energy electrons and holes; p,(¢,) and p,(¢,) 
are the concentrations of electrons and holes fixed in 
traps with a depth of ¢, (¢,); 4, and ys, (4, and pp) are 
the mobilities of the high-energy (low-energy) electrons 
and holes, respectively; k,, &,, and &, are electron-hole 
recombination rate constants; 6, and (4, are the rate 
constants of capture of electrons and holes in the traps; 
k,' (ky) are rate constants of the release of electrons 
(holes) from traps with a depth of ¢, (€,); g is the density 
of the electrical charge; 7 is the temperature; x is the 
dielectric constant; a and § are quantities that are the 
inverse of the lifetime of the high-energy electrons and 
holes'”; M,(¢,) and M,(¢,) are the concentrations of traps 
for electrons and holes with a depth of «, (¢,); and K,, and 
K, are the number of electron and hole energy levels in 
the immobility band that are being examined. 
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Relaxation of the space charge is accompanied by the 
release of heat (Q,): 


Q=aBV, cme> uA, 


where @ is radiation conductance and ju, and n, are the 
mobility and concentration of the i-th type of charge 
Carrier. 


The function G(h,t) is related to the function Q,(h,1) by 
the following relationship: 


G(h, 1 = 2% 9/ Qe 


where Q., is the energy expended to create one charge 
carrier pair. The average amount of energy expended 
to create one electron-hole pair is approximately equal 
to 1.5¢,, with the overwhelming majority of irradiation- 
generated carriers being high-energy carriers.’ They are 
thermalized within a short time. Not all of the carriers 
initially formed are capable of taking part in charge 
transfer, however. A significant portion lose their kinetic 
energy before they have time to move beyond the 
confines of their partner's sphere of Coulomb attraction 
and recombine.® As a result, the effective value of Q_. is 
much larger than 1.5¢,. When an external electrical field 
is applied and the temperature is increased, the yield of 
free carriers should increase.*° It this study it has been 
assumed that Q., = 100 eV. ; 


It has been assumed that the traps for electrons (holes) 
are structured and specialized and that their concentra- 
tions, i.c., M,(¢,) (M,(€,)), are exponentially dependent 
on their depth «, (¢,)'*: 


M, («) = M, exp (- 8/&7)), 
M, (8) = M, exp (- e/kT,), 


where T, and T, are characteristic constants for the given 
material; M, and M, are determined by the normalization 
conditions; UM,(¢,) = C, and LM fe) = C,; C, and C, 
are the total concentrations of traps for electrons (holes); 
and & is Boltzmann's constant. 


The dependence of the coefficients in equation system 
(2) on temperature and the electric field and the rela- 
tionships between them that have been assumed in this 
study are reflected in Table 1. 
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Table 1. 
m0) Seemoeecn (2) marge 
Ae Ky Coasts re 
or By ar {1, 10) 
by hy hy hy Hn hy ign by eH (1, 8, 11) 
(4)@ yeovers ftp >> p,, (3) 
Bar Bp By “HyT. Bp “Hg? G) (1, 8, 11) 
Mp i, | senp(—(4, ~ anf)/k7), om e— ¢, — exyhom soayemn (12) 


i = coefficient; 2 = relationship; 3 = reference no.; 4 = with 
Consideration for; 5 = half-width; and 6 = trap depth. 


Equation system (2) must be supplemented with two 
types of boundary conditions. 

The first is a condition regarding the electrical field 
that is dictated by the problem's geometry (the presence 
and location of grounded electrons," irradiation of the 
study object in a vacuum or in a medium with high 
conduction, etc.). In our case, this condition has the 
following form!: 


a 
f B(x) dx = 0 


(H being the distance between electrodes). 


The second is the condition of exchange of charge car- 
riers between the dielectric and ing medium 
(metal electrode, vacuum, plasma, etc.).'*'’ The condi- 
tion of blocking electrodes' is used in this study. The 
condition is as follows: The charge carriers may leave 
the study object unimpeded; however, they cannot be 
injected inward from the electrons. 

The mathematical model described above was used to 
calculate the effect that a monoenergetic relativistic 
electron beam with parameters typical for modern small 
accelerators (an electron energy in the beam of 0.5 
MeV, an average current density of 0.5 kA/cm’, and a 
sinusoidal pulse with a duration of 50 ns) has on a 2.5- 
mm-thick polystyrene wafer located inside a grounded 
metal cladding. 
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Figure 2. Distributions of electron charge density 
throughout the thickness of the study object at the 
following moments in time: 1, 5 ns; 2, 10 ns; and 3, 
1,500 ns. (Values along the x-axis represent é in 
cm, and values along the y-axis represent @ in 
coulombe/cm’.) 


The calculations were performed twice: with (case A) 
and without (case B) consideration for the processes of 
space charge accumulation and relaxation. 
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Figure 3. Distributions of the electrical field 
throughout the thickness of the study object at the 
following moments in time: 1, 5 ns; 2, 10 ns; and 3, 
1,500 ns. (Values along the x-axis represent h in cm, 
and values along the y-axis represent E in MV/cm.) 


Figures 2 through 8 present the results of the calcu- 
lations. The beginning of the irradiation process was 
marked by rapid accumulation of a space charge in the 
study object (Figure 2), accompanied by an increase in 
the internal electrical field (Figure 3). After approxi- 
mately 10 ns, when the dielectric’s conduction had in- 
creased significantly (under the stated conditions, from 
10° to 10* ohm/cm), the said process nearly halted, af- 
ter which the distributions of the charge and the electric 
field in the study object remained virtually unchanged. 
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Figure 4. Distributions of specific internal energy 
throughout the thickness of the study object (curves 
1 and 2 illustrate case A, whereas curves 3 and 4 
illustrate case B) at the following moments in time: 
1 and 3, 25 ns; 2 and 4, S@ ns. (Values along the 
x-axis represent h in cm; values along the y-axis 
represent ¢ in J/g.) 
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From that moment, virtually the entire charge injected 
by the beam into the study object flowed out through its 
anterior surface. This means that a conduction current 
with a density on the order of the current density was 
induced in the study object and accompanied by strong 
ohmic heating of the object. We will note that under 
similar conditions in metals, ohmic heating may often 
be ignored. 
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Figure 5. Time dependences of the velocity of the 
study object’s free surface: 1, case A; 2, case B. 
(Values along the x-axis represent ¢ in ns, and 
values along the y-axis represent U x 10° km/s.) 


During the irradiation process, between 20 and 30 
percent of the beam’s total energy is released in the 
form of Joulian heat. The mechanism of redistribution 
of this energy is as follows: Upon reaching the target, 
a fast electron expends approximately 75 percent of its 
energy on ionization (excitation) of the molecules of the 
matter and approximately 25 percent on completing the 
work involved in overcoming the repulsive force of the 
negative charge that has accumulated in the study object. 
It is obvious that the potential energy of the charge in 
the target increases by an amount equal to that work 
and that the said energy is, in turn, released virtually 
immediately in ohmic heating of the study object by 
conduction currents. We will note that the total electrical 
energy in the study object equals 


a 


where integration is performed with respect to the study 
object's entire volume and that at any given moment in 
time, it does not exceed | percent of the beam’s total 
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Figure 6. Distributions of the concentration of free 
low-energy electrons throughout the thickness of 
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1, 5 ms; 2, 10 ms; 3, 25 ns; and 4, 51 ns. (Values 
along the x-axis represent h in cm, and values along 
the y-axis represent NV, x 10'°/cm’.) 


The occurrence in the study object of an electrical field 
and a conduction current induced by the electrical field 
has a significant effect on the spatial distribution of 
the energy contribution (see Figure 4): The area of 
energy release narrows (the electrical field of the space 
charge exerts a stopping action on the electrons of the 
beam), the profile of the energy release curve in the 
near-surface area of the study object changes shape 
(including shape changes caused by ohmic heating), and 
the energy density in the study object increase. This 
results in an increase in absolute voltage in the study 
object and the velocity of its free surface. In addition, 
narrowing of the region of the energy contribution (by 
approximately 30 percent) leads to a marked slowing 
of the movement of the compression wave out onto the 
rear surface of the target (Figure 5). 
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Figure 7. Distributions of the concentration of 
electrons fixed in traps throughout the thickness of 


the study object at the following moments in time: 
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1, 5 ns; 2, 10 ns; 3, 25 ns; 4, 51 ns; and 5, 1,500 ns. 
(Values along the x-axis represent 4 in cm, and 
values along the y-axis represent Nx 10'*=’.) 


The dynamics of accumulation and relaxation of the 
space charge is especially interesting. The charge is 
implanted by the beam primarily in the middle part of 
the target; consequently, under the effect of the electrical 
field, it moves toward the anterior surface. From there, 
it flows down onto the grounded cladding and into 
the depths of the study object where it is accumulated 
(see Figure 2) because the sharp decrease in radiation 
conductance in the object's depths, which exceeds the 
pathlength of beam’s electrons, results in the virtually 
total absence of current flowing to the anterior surface 
of the study object. Despite the fact that a large part 
of the charge injected into the dielectric flows through 
the anterior surface, part of the charge accumulates in 
the sample, which in turn leads to the development 
of a strong electrical field that is capable of causing 
breakdown of the material (see Figure 3). 
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Figure 8. Distributions of the concentration 
throughout the thickness of the study object at the 
following moments in time: 1, 5 ns; 2, 10 ns; 3, 25 
ns; 4, 45 ns; and 5, 51 ns. (Values along the x-axis 

represent h in cm, and values along the y-axis 
represent n x 10'°/cm’.) 


The kinetics of the charge carriers is characterized by 
the following features (see Figures 6 through 8): During 
the effect of the relativistic electron beam, the free low- 
energy Cafriers in the study object constitute approxi- 
mately 93 percent of the total numbe. of charge carriers. 
Only approximately 6 percent of the carriers are fixed 
on traps, and the high-energy carriers (with the given 
parameters of the relativistic electron beam) amount to 
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about | percent. After the irradiation has ended, the pic- 
ture changes: The high-energy carriers disappear alto- 
gether, and the concentration of free low-energy carriers 
drops sharply; however, a significant number of charge 
carriers remain stabilized on traps in the study object. 
Because recombination between the electrons and holes 
fixed on traps is impeded (it is not taken into consid- 
eration within the framework of the given model) and 
the number of free carriers is insignificant, the said im- 
balanced excess concentration of charge carriers in the 
study object may persist for tens of hours. In particu- 
lar, thermal ejection of the carriers from the traps in the 
immobility band results in slow discharge of the study 
object (postradiation conductance) and anomalous heat- 
stimulated conduction.'! 


All of the aforesaid applies to the case where no 
electrical breakdown occurs in the study object during 
the irradiation. In the opposite event, the pattern of 
the interaction of a relativistic electron beam with a 
dielectric would be fundamentally altered. As a result 
of breakdown, the study object is discharged and its 
internal electrical field drops virtually to zero. It is 
obvious that no effects related to ohmic heating and 
contraction of the energy release band may be present. 
The energy that is released in the breakdown channel 
will equal the electrical energy «, accumulated ia the 
sample at the moment when the discharge begins, 
i.e., as was indicated above, it constitutes less than | 
percent of the beam’s total energy. Consequently, the 
disturbances that the discharge causes in the total picture 
of the hydrodynamic flow of the medium irradiated by 
a relativistic electron beam should not be very great. 
Because the aforesaid energy will be released in a 
small space (on the order of the total volume of the 
breakdown channels), however, local disturbances, such 
as a deterioration of the plane geometry of the shock 
wave or melting through of the study object along the 
branches of the breakdown tree (which in this case 
should be initiated in the depths of the sample and 
extend to its anterior surface), etc., may take place. 
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(FBIS Translated Text) Numerical calculations are used 
to study the effect of backpressure on drag coefficient 
when a sphere penetrates a liquid in the presence of 
backpressure. It is demonstrated that backpressure with 
cavitation numbers of o > 1.5 results in a decrease in 
drag coefficient. Analytic solutions are used to demon- 
strate that taking the decrease in drag coefficient into 
account affects estimates of the degree of disordering of 
the material that is required to explain ultradeep pene- 
tration. 


Introduction 


Different mechanisms that could explain the phenome- 
non of ultradeep penetration have been widely discussed 
in the literature.'* The point of departure for our study 
was a publication by Altshuler et al.‘ that hypothesized 
that two factors are mandatory for development of the 
phenomenon of ultradeep penetration: brief disordering 
of tne material of the barrier on account of uneven heat- 
ing under conditions of ultrarapid deformation and the 
development of a mode of forced laminar circulation in 
the disordered medium ("quasi-liquid") under the effect 
of external pressures. 


To estimate the effect of these factors, a numerical study 
of the motion of a body in a medium with backpressure 
was performed,’ and analytic estimates of the degree of 
disordering were performed.° 


This article combines the two aforesaid studies, , and 
presents additional material illustrating the process of 
development of forced quasi-laminar circulation and a 
decrease in the drag coefficient. 


1. Numerical Investigation of the Motion of a 
Sphere in a Medium With Backpressure 


1) The problem of the motion of a nondeformed sphere 
R, = 1 in a condensed medium was examined in a hy- 
drodynamic approximation in a two-dimensional for- 
mulation by using a method that has been described 
elsewhere’* and that uses finite-difference approxima- 
tions of mass, motion, and energy conservation laws. 
The medium's equation of state was assumed to have a 
Mie-Grunhausen form’ with the following values: p, = 
7.8 gicm’, C, = 5 km/s, n = 3, and [' = 2.548. The vari- 
ous values of 


values. The motion of the entire medium surrounding 
the sphere was specified by the velocity u, = 1 km/s, 


which corresponds to the motion of a sphere with a ve- 
locity of -a, (with a Mach number (M) equal to 0.2). 


FBIS-UST-95-052 
19 December 1995 


The number of points in the main series of calculations 
was as follows: 40 in a radial direction and 34 at an 
angle. The points in a radiation direction were selected 
$0 as to be denser toward the sphere’s boundary with a 
denominator of g = 0.94. The size of the point adjacent 
to the sphere equaled 0.05. 


The force of the sphere’s drag F in the liquid was 
calculated by integration of the quantity pdS along the 
sphere’s surface: 


f= Sods, = -c, 2S ari, 


where R, is the sphere’s radius. The dependence of the 
drag coefficient (C,] on the magnitude of backpressure 
(p,,) and cavitation number [0] was examined in the 
calculations. 


The literature'*'' contains data regarding the effect of 
relatively small cavitation numbers, i.e., 0 = 2p_/p_u,? 
(where p, and p, are the pressure and density in the 
medium and u is the velocity of the oncoming flow), on 
the magnitude of the drag coefficient [C,}. In this study, 
the effect of pressure is investigated numerically within 
a wider range. 


2) Figures 1 and 2 and Table | present the main results 
of the calculations. 
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Figure 1. Dependence of drag coefficient for 


different cavitation numbers: 1, o = 0; 2, o = 0.5; 


Figure | (curve |) presents the function C,(s) for a 
cavitation number of zero. After the maximum value 
C, nex 88 been reached, the value of C, decreases and 
reaches some steady-state value C, . once a steady state 


3, o = 16.7. 
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has been established. The value C, = 0.37, which was 
obtained in calculations for a steady state, is in rather 
good agreement with data of approximate theory,'®'! 
according to which C, = 0.34. Table | presents the main 
calculations and selected results. 
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Figure 2. Dependence of the drag coefficient in a 
steady flow on the cavitation number. 





Figure 2 presents the calculated dependence of C, , 
on the cavitation number. In cases of large cavitation 
numbers (0 < 0.8), the calculated data are consistent 
with the published'*'' linear dependence of C, on the 
cavitation number: C (7) = C.(0\1 + @). 





lanes 


Table 1. Drag Coefficient for Different Cavitation Numbers 









































Initial V abwes Drag Coefiictent 

Ne. o P,,» GPa P» @/cm’ = Cy 

! 0.0 00 78 6.82 0.37 
2 0.166 0.65 7.826 7.0 0.46 
3 0.495 1.95 7.877 7.37 0.58 
4 0.820 3.25 7.928 7.74 0.66 
$ 0.981 3.9 7.953 7.92 0.50 
6 1.220 4.875 7.990 8.18 0.31 
? 1.458 5.85 8.027 8.44 0.1 

” 1.926 78 8.100 8.95 0.09 
9 2.385 9.75 8.175 9.18 0.09 
10 10.37 48.75 9.400 14.797 0.09 
i" 18.41 97.5 10.59 16.68 0.09 
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The flow velocity fields and cavern boundaries pre- In the case of zero cavitation (Figure 3), a separation 
sented in Figures 3 through 5 provide a graphic illustra- point (¢ = 58°) and cavern on the rear side of the body 
tion of the nature of steady-state circulation with cav- are evident. 

itation numbers of a = 0, 0.82, and 16.7. The cavern 

boundary was plotted as the isobar corresponding to a 

pressure of 0.05 GPa. 








Figure 3. Velocity field and position of the cavern in a steady flow with zero backpressure and a cavitation 
number of o = 0. 


When the cavitation number is increased (Figure 4), changed. With cavitation numbers of o = 0.5 and 0.82, 
the separation point is shifted, the velocity field is the separation points correspond to angles of ¢ = 69° 
restructured, and the configuration of the cavern is and 74°, respectively. 
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Figure 4. Velocity field and position of the cavern in a steady flow with a cavitation number of o = 0.82. 


The velocity field corresponding to a cavitation number _-velocity field of a laminar flow, which corresponds to a 
of o = 16.7 (Figure 5) is virtually identical to the decrease in the drag coefficient. No cavern is present. 














: J + 
Figure 5. Velocity field in a steady flow with a cavitation number of o = 16.7; no cavern is present. 


The increase in backpressure occurring with cavitation  Backpressure complicates the flow of the medium in 
sumbers of o > 0.8 thus results in a decrease in C,.. It a nonsteady state. At specified moments in time, the 
is evident from Table | and Figure 2 that with cavitation collapse of the cavern causes the pressure on the half- 
sumbers of 7 > 1.5, C_ = 0.09, which is one-fourth the sphere opposite the oncoming flow to become greater 
value for the case where the cavitation number equals _— than that on the half-sphere turned toward the oncoming 
zero. flow. The drag coefficient C_ becomes negative. This is 
seen in Figure 1, where curves |, 2, and 3 represent the 
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function C(t) for cavitation numbers of a = 0, 0.5, and 
16.7, respectively. 
Table 2 illustrates the development of a negative drag 


coefficient. Table 2 presents the pressures at the polar 
points of the sphere (¢ = 0 and 180°) and at the equator 
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(@ = 90°) at different moments in time (7 = 0.39, 
0.68, and 1.8) for a cavitation number of 16.7, which 
corresponds to a backpressure of 65 GPa. The drag 
coefficients (C_) corresponding to those times may be 
obtained by examining Figure |. 





Table 2. Pressure (in GPa) at Points on the Sphere 

















¢ 
T 0° 90° 180 
0.39 67.5 $7.2 73.8 
0.68 68.7 59.5 69.7 
1.8 70 59.5 67.4 














It is especially important to emphasize that the negative 
drag coefficient occurring with large cavitation numbers 
in a nonsteady state cannot explain ultradeep penetration 
resulting from the brief action of the force accelerating 
the body. The body covers a distance on the order 
of a radius during the existence of the negative drag 
coefficient. The calculated decrease in drag coefficient 
in the steady-state stage of penetration from C, = 0.37 
when o = 0 to C, = 0.09 with cavitation numbers 
of a > 1.5 is important for the phenomenon under 
consideration. 


2. Estimation of the Required Degree of 
Disordering of the Medium 


1) The case of an elastoplastic medium. The force acting 


on a body moving in an elastoplastic medium may be 
represented in the following form'?-"’: 


d 
FamG=-c,%s- Hy, (4) 


where H is the material's dynamic hardness, which 
depends on the yield strength Y and equals H = 3Y, S is 
the cross-sectional area of the middle, and the remaining 
designations are generally accepted. 


By integrating (1), we obtain 


u (0) = (b/a)* tg (9, - ¢ (ed), 


em arcs (F) "uy oe SE, 6 = SS. 


It follows from (2) that when f° = y/(ab)'”, the body 
stops. The penetration depth L(t) is determined by the 
relationship 


£ (0) = Lin toon {¥% ~ #(08)1 - +10 100s yl. 


At the moment the body stops, we have 


b= bin too wi"! = oe la (1 + SPB). (3) 


For the case of a sphere with a density of p, and 
diameter of d, we obtain 


26 
Laie 


in (1+ 4), (4) 


where L* = Lid and Y° = 2Y/pu,? are the reduced 
depth and yield strength [Note: here and henceforth, the 
asterisk following the variables L and Y is used in place 
of a tilde) and 6 = p,/p is the ratio of the densities of 
the body and medium. 
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It follows from relationship (4) that when Y » C., the 
penetration depth is virtually independent of C,. In that 
case, L* << | and strength properties are decisive. 
From relationship (4) we determine Y° as a function of 
L’: 


37 


with backpressure will have on the magnitude of the 
yield strength with different penetration depths. 


Table 3 presents the calculation results. It was assumed 
that 6 = 1, p = 7.8 g/cm’, and u, = | km/s. It is 
evident from an examination of the values presented in 
Table 3 that a decrease in C, has a significant effect on 
the permissible yield strength values in cases of large 























Fa S [exp (3:2) -17". (5) penetration depths (L* 2 50). 
Relationship (5) makes it possible to analyze the effect 
that a change in the drag coefficient (C,) in connection 
Table 3. 
T Y (kbar) when C, = 0.37 (p = 0) Y (kbar) when C, = 0.09 (0 > 1.5) 
1 8.1 
2 3.8 
10 0.019 0.41 
50 4.27x 10°" 1.37 10° 
100 3.79 x 10% 16x 10° 











2) The case of a viscous medium. The force acting on 


a sphere moving in a viscous liquid may be represented 
in the following form'*: 


mdu 
yen @ ~C,2E 5 - any du, (6) 
where ys is the dynamic viscosity coefficient. 


Integration of equation (6) results in the following 
expression for the penetration rate, i.e., u(t) and depth, 


ie., L(t): 
a 
80- sega a 
L(g=tm[i+Sa-«9). 


where g = 3nyd/m. It follows from the preceding 
relationship that the maximum penetration depth equals 


be fia [i+]. 


We thus obtain the following expression for the reduced 
penetration depth: 


LaZdin (i+ spre), an 


where Re = udp/j is the Reynolds number. 

When Re << 2 (high viscosity), L° is weakly dependent 
on C,. Everything is determined by the viscosity rather 
than the velocity head. 

By using relationship (7), we represent the dynamic 
viscosity coefficient (4) as a function of penetration 
depth: 


ne Se [om (252) -1)". 
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Just as in the case of an clastoplastic medium, we will 
trace the effect of a change in C, on the magnitude of the 
dynamic viscosity coefficient, which must be ascribed 
to the medium at different penetration depths. Table 4 
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presents the calculation results. When the calculations 
were made, it was assumed that 6 = 1, p = 7.8 g/cm’, 
u, = | km/s, and d = 20 zm = 2 x 10° cm. 























_ _— 
Table 4. 
L’ js (pt) whea C, = 0.37 (p = 0) js (pa) when C, = 0.09 (0 > 1.5) 
1 0.75 x 10 0.84 x 10° 
2 0.32% 10 04x10 
10 6 
50 2.26 x 10° 2.07 x 10° 
100 2.13 x 10°"! 6.86 x 10° 











Note: pz = pieze; L” = the variable L with a tilde, which occurs in displayed equations throughout this article. 





It is evident from the data in Table 4 that just as in 
the case in an elastoplastic medium, a change in C, 
significantly alters the permissible values of the dynamic 
viscosity coefficient, which must be ascribed io the 
medium at large penetration depths (2 50). 


The numerical computations that were performed con- 
firmed the published* hypothesis regarding the possibil- 
ity of development of forced quasi-laminar circulation 
and a decrease in drag coefficient during motion in a 
medium with backpressure. It has been demonstrated 


that this effect occurs at pressures corresponding to cav- 
itation numbers of o > 1.5. 


Analytic solutions have helped demonstrate that taking 
the drop in the drag coefficient to one-fourth its value 
(from 0.37 to 0.09) into consideration has a significant 
effect on estimates of the degree of disordering of the 
material that is required to explain ultradeep penetration. 
The authors thank L.V. Altshuler, at whose suggestion 
and with whose support this study was undertaken. 
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